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ABSTRACT 


In  the  linear  theory  of  elasticity,  Saint- Venant’s  principle  is  used  to  justify  the  ne¬ 
glect  of  edge  effects  when  determining  stresses  in  a  body.  For  isotropic  materials,  the 
validity  of  this  is  well  established.  However  for  anisotropic  materials,  experimental 
results  have  shown  that  edge  effects  may  persist  much  farther  into  the  material  than 
for  isotropic  materials  and  as  a  result  cannot  be  neglected.  This  research  examines 
the  effect  of  material  anisotropy  on  the  exponential  decay  rate  for  stresses  in  a  semi¬ 
infinite  elastic  strip.  A  linear  elastic  semi-infinite  strip  in  a  state  of  plane  stress/strain 
subject  to  a  self-equilibrated  end  load  is  considered  first  for  a  specially  orthotropic 
material  and  then  for  the  general  anisotropic  material.  The  problem  is  governed  by 
a  fourth-order  elliptic  partial  differential  equation  with  constant  coefficients.  Con¬ 
servation  properties  of  the  solution  axe  derived  to  help  to  determine  a  stress  decay 
rate  estimate.  Energy  methods  are  then  used  to  establish  lower  bounds  on  the  actual 
stress  decay  rate.  Both  analytic  and  numerical  estimates  are  obtained  in  terms  of 
the  elastic  constants  of  the  material  and  results  are  shown  for  several  contemporary 
engineering  materials.  When  compared  with  the  exact  stress  decay  rate  computed 
numerically  from  the  eigenvalues  of  a  fourth-order  ordinary  differential  equation,  the 
results  in  some  cases  show  a  high  degree  of  accuracy.  Results  of  the  type  obtained 
here  have  several  important  practical  applications.  For  example,  the  results  provide 
physical  insight  into  the  mechanical  testing  of  anisotropic  and  laminated  composite 
structures  and  are  useful  in  assessing  the  influence  of  fasteners,  joints,  etc.  on  the 
behavior  of  composite  structures. 
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Chapter  1 


INTRODUCTION 


When  considering  plane  deformations  in  the  context  of  linear  elasticity  theory,  Saint- 
Venant’s  principle  is  often  used  to  justify  approximations  that  neglect  edge  eflFects 
(see  [1],  [2]  for  recent  reviews).  As  a  result,  this  principle  is  also  the  basis  for  many 
approximate  “strength  of  materials”  formulae  that  are  widely  used  in  engineering 
practice.  Experimental  evidence  has  supported  the  use  of  Saint- Venant’s  principle  in 
this  manner.  For  instance,  it  is  commonly  known  that  for  structures  made  of  homoge¬ 
neous  isotropic  materials,  stress  effects  due  to  self-equilibrated  edge  loads  decay  with 
distance  from  the  loaded  ends  over  relatively  short  lengths,  e.g.  one  strip  width  for  a 
semi-infinite  strip.  For  structures  made  of  homogeneous  anisotropic  materials,  how¬ 
ever,  it  has  been  shown  that  these  local  stress  effects  may  persist  well  into  the  interior 
of  the  body  and  thus  cannot  be  neglected  when  making  approximations  (see  [1]  -[10]). 
It  is  the  goal  of  this  dissertation  to  investigate  the  behavior  of  local  stress  effects  for 
anisotropic  materials  with  a  particular  emphasis  on  determining  the  influence  of  the 
elastic  constants  of  the  material  on  the  rate  of  stress  decay. 
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The  paxticulax  problem  that  we  consider  is  that  of  a  linearly  elastic  semi-infinite 
strip  in  a  state  of  plane  stress/strain  subject  to  a  self-equilibrated  end  load.  For  strips 
made  of  specially  orthotropic  materials,  this  stress  decay  problem  may  be  formulated 
in  terms  of  a  single  dimensionless  material  parameter.  For  strips  made  of  general 
anisotropic  materials,  however,  the  problem  is  more  complex  and  involves  three  di¬ 
mensionless  material  parameters.  In  both  cases,  the  elasticity  problem  is  governed  by 
a  fourth-order  elliptic  partial  differential  equation  with  constant  coefficients,  whose 
exact  analytical  solutions  are  difficult  to  obtain.  In  studying  the  effects  of  material 
anisotropy  on  the  exponential  decay  for  stresses,  we  will  develop  methods  of  exam¬ 
ining  the  relevant  elasticity  problems  without  having  to  find  exact  solutions.  Energy 
methods  and  estimating  techniques  for  partial  differential  equations  will  be  used  to 
obtain  bounds  on  the  solution  and  on  the  stresses  (which  are  second  derivatives  of 
the  solution)  with  an  emphasis  on  obtaining  bounds  for  their  decay  rates.  In  addition 
to  the  specific  problem  of  plane  elasticity  for  a  semi-infinite  strip,  these  methods  can 
also  be  applied  to  other  more  general  situations,  i.e.  problems  with  different  geome¬ 
tries,  problems  with  different  governing  equations,  nonlinear  problems,  etc.  (see  e.g. 
[1],  [2]).  Through  the  use  of  such  energy  methods  and  estimating  techniques,  both 
analytic  and  numerical  decay  rate  estimates  axe  obtained  in  terms  of  the  elastic  con¬ 
stants  of  the  material.  For  a  set  of  specific  materials,  these  estimates  are  compared 
with  the  exact  decay  rate  computed  numerically  from  the  eigenvalues  of  a  fourth- 
order  ordinary  differential  equation.  The  results  in  some  cases  show  a  high  degree  of 
accuracy. 

Much  work  has  already  been  done  with  isotropic  materials  in  a  state  of  plane 
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stress/strain  for  which  the  governing  equation  is  the  bihaxmonic  equation  (see  [3], 
[11]  -[15]).  Various  energy  methods  and  differential  inequality  techniques  have  been 
developed  to  obtain  lower  bounds  for  the  decay  rate  that  accurately  estimate  the 
exact  decay  rate  of  the  solution.  Fewer  results  eire  known  for  anisotropic  materials  in  a 
state  of  plane  stress/strain  where  the  governing  equation  is  a  generalized  fourth-order 
elliptic  partial  differential  equation  with  constant  coefficients.  Much  of  the  analysis 
that  has  already  been  done  for  anisotropic  materials  has  been  concerned  with  exact 
numerical  solutions  to  the  problem  rather  than  analytic  estimates.  This  dissertation 
will  extend  some  of  the  known  techniques  used  for  the  biharmonic  equation  and  apply 
them  to  the  generalized  fourth-order  problem.  In  the  process,  new  techniques  are  also 
generated.  Ultimately,  an  explicit  “formula”  for  the  decay  rate  of  stresses  is  sought 
which,  although  it  will  be  an  estimate  (in  fact,  a  lower  bound  for  the  exact  decay 
rate),  will  reveal  how  the  elastic  constants  are  involved. 

This  research  has  several  applications  in  the  field  of  structural  analysis  and  design. 
For  example,  in  the  mechanical  testing  of  anisotropic  and  composite  materials,  regions 
of  uniform  stress  and  strain  fields  are  necessary  for  accurate  measurements  of  material 
properties.  Thus  it  is  crucial  to  have  an  understanding  of  the  local  stress  behavior 
and  edge  effects  that  may  be  induced  through  the  testing  procedure.  Another  area  of 
application  is  in  assessing  the  influence  of  fasteners,  joints,  rivets,  etc.  in  composite 
materials.  Here,  an  understanding  of  local  stress  effects  and  how  fax  into  the  material 
they  penetrate  is  of  critical  importance  to  the  designer.  Results  of  the  type  obtained 
in  this  research  could  also  allow  for  “tailoring”  a  material  with  specific  properties  to 


ensure  that  local  stresses  attenuate  at  a  desired  rate. 
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The  organization  of  this  thesis  is  as  follows.  In  Chapter  2,  we  formulate  the 
anisotropic  strip  problem  and  its  nondimensional  form.  Several  conservation  prop¬ 
erties  are  derived  for  later  use.  The  strain-energy  for  the  problem  is  examined,  the 
positive-definiteness  of  which  leads  to  restrictions  on  the  material  constants.  Energy 
methods  that  involve  estimating  techniques  using  differential  inequalities  are  outlined 
and  several  energy  norms  are  presented.  In  Chapter  3,  we  examine  these  energy  meth¬ 
ods  for  a  special  class  of  materials,  namely  orthotropic  materials  and  the  problem  is 
reduced  to  the  orthotropic  strip  problem.  This  serves  as  a  simpler  setting  within  which 
to  illustrate  the  energy  techniques  of  this  dissertation.  Several  decay  rate  estimates 
are  developed  and  then  compared  with  exact  solutions.  Where  possible,  explicit  for¬ 
mulae  for  the  estimated  decay  rates  are  obtained  and  the  results  are  compared  with 
exact  decay  rates  (computed  from  roots  of  a  transcendental  equation)  for  a  variety  of 
materials  used  in  composites  technology.  Asymptotic  results  for  strongly  orthotropic 
materials  are  also  developed.  In  Chapter  4,  we  return  to  the  anisotropic  problem 
and  extend  the  techniques  from  Chapter  3  to  cover  more  general  materials.  The 
complexity  of  the  analysis  is  greatly  increased.  Several  estimates  are  developed  and 
numerical  results  are  compared  to  exact  solutions.  The  accuracy  of  these  results  and 
their  implications  for  future  estimating  techniques  are  discussed.  Finally,  in  Chapter 
5  we  summarize  the  important  results  obtained  in  this  dissertation  and  indicate  some 


directions  for  future  research. 


Chapter  2 


ANISOTROPIC  ELASTIC 
STRIP  PROBLEM 


In  this  chapter,  the  anisotropic  elcistic  strip  problem  is  formulated  and  preliminary 
material  needed  for  the  estimation  techniques  used  in  later  chapters  is  presented.  The 
chapter  begins  with  a  formulation  of  the  physical  problem,  which  is  then  nondimen- 
sionalized.  Next,  several  conservation  properties  are  examined  as  well  as  positive¬ 
definiteness  of  the  strain-energy  density.  The  chapter  ends  with  an  outline  of  the 
energy  arguments  that  form  the  basis  for  subsequent  estimation  techniques. 

2.1  Problem  Formulation 

Consider  an  homogeneous  anisotropic  linearly  elastic  body  that  occupies  the  following 
region  described  in  the  xi,  X2,  X3  Cartesian  coordinate  system  by 

xi  >0,  -H  <  X2  <  H,  -T  <  X3  <  T,  (2.1) 
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where  2H  is  the  height  of  the  body  and  2T  is  the  thickness  that  may  be  taken  to 
be  infinite  or  infinitesimal.  For  a  material  with  the  xi  —  X2  plane  a  plane  of  elastic 
symmetry,  one  can  define  a  state  of  plane  deformation  such  that  the  infinitely  thick 
body  is  in  a  state  of  plane  strain  or  such  that  the  infinitesimally  thin  body  is  in  a 
state  of  plane  stress.  All  surfaces  of  the  body  will  be  traction  free  except  for  the 
surface  xi  =  0,  where  a  prescribed  self-equilibrated  traction  is  applied.  For  this  type 
of  loading,  it  is  assumed  that  the  stresses  decay  to  zero  as  xi  — >  oo.  The  two  plane 
strain  and  plane  stress  elastostatic  problems  are  identical  mathematically  and  differ 
only  in  the  numerical  value  of  the  elastic  constants  used  for  each  problem.  Both  result 
in  the  consideration  of  an  anisotropic  elastic  semi-infinite  strip  occupying  the  region 
3i  in  the  Xi  —  X2  plane  defined  by 

:  Si  >  0,  -H  <  X2  <  H,  (2.2) 

with  traction-free  lateral  surfaces  ,  a  prescribed  self-equilibrated  edge  load  on  the  end 
Si  =  0  and  stresses  that  decay  to  zero  as  Si  — »•  oo  (see  Figure  2.1). 


Figure  2.1:  Semi-infinite  strip 

The  in-plane  components  of  the  strain  tensor  en,  ei2,  622  axe  related  to  the  in- 
plane  components  of  the  stress  tensor  th  ,  T12,  T22  by  the  following  reduced  constitutive 
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equations  (see  e.g.  [4]), 


fill  =  i^llTll  +  ^12^22  +  fil6Tl2i 

622  =  021^11  +  /I22T22  +  ^26^12, 

2ei2  =  ^61 'Til  +  ^62^22  +  ^66^12) 


or  in  matrix  notation 


en 

fin  fin  file 

Til 

622 

fi21  fi22  fi26 

T22 

26i2 

fiei  fie2  fiee 

T12 

(2.3) 


(2.4) 


The  constitutive  equations  are  obtained  from  the  general  3-D  constitutive  equations 
by  invoking  the  a,ssumptions  of  plane  strain  or  plane  stress  to  obtain  a  reduced  2-D 
form.  The  constants  =  ^gp  (p,g  =  1,2, 6)  can  be  written  in  terms  of  the  elastic 
compliances  Opg.  In  the  case  of  plane  stress,  the  fipg’s  are  the  elastic  compliances 
themselves,  while  for  plane  strain  they  are  related  by 

^pg  =  apg  -  (033  0).  (2.5) 

033 

For  materials  that  are  transversely  isotropic  about  the  Xi  axis  (or  for  specially  or¬ 
thotropic  materials  whose  fibers  run  in  the  x\  direction),  considerable  simplification 
results  from  the  fact  that 

022  =  ®33j  ®12  =  ®13,  <*16  =  <*26  =  <*36  =  0-  (2*®) 


As  a  result,  /3i6  =  ^26  =  0  and  the  elastic  constants  ^P9  may  be  expressed  in  terms 
of  the  usual  engineering  constants  v,  E  and  G  which  denote  Poisson’s  ratio,  Young’s 
modulus  and  shear  modulus  respectively,  so  that  one  obtains 
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plane  stress  :  = 


^22  = 


J_ 

El' 

J_ 

Et' 


= 


VLT 

'  El 


Glt' 


^16  =  ^26  =  0 


(2.7) 


plane  strain  :  jin  =  "^(l  —  v\jEtI El)i 


^22  = 


El' 

—(1  —  I'xr), 
£/T 


Pl2  = 

1 

Glt' 


^LT 

El 


(1  +  vtt) 


^16  =  ^26  =  0.  (2.8) 


In  these  expressions,  L  denotes  the  longitudinal  direction  parallel  to  the  xi  axis  and  T 
denotes  the  transverse  direction.  Further  simplification  arises  in  the  case  of  isotropic 
materials  such  that 


plane  stress 


1  V 

Pn  —  P22  =  Pl2  =  — ^ 

^66  =  -Q,  ^16  =  /?26  =  0 


(2.9) 


1  2  ^ 

plane  strain  :  /3n  =  1^22  —  t;(1  -  i'  ),  ^12  =  -^(1  +  v) 

h/  h/ 

^66  =  -Q,  /3i6  =  ^26  =  0.  (2.10) 

Returning  to  the  anisotropic  analysis,  an  Airy  stress  function  ^(xi,X2)  is  defined 
in  terms  of  which  the  Cartesian  components  of  the  stress  tensor  r  are  given  by 

Til  =  ^,22j  T22  =  T\2  =  —<i>,\2  J  (2-11) 

where  a  comma  denotes  partial  differentiation.  With  the  Airy  stress  function  defined 
in  (2.11),  the  equations  of  equilibrium  (with  zero  body  force) 
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are  identically  satisfied.  The  compatibility  equations  then  yield  the  following  govern¬ 
ing  differential  equation  on  (see  e.g.  [4],  [16]), 

^22^,1111  —  2^26^,1112  +  (2^12  +  ^66)^,1122 

—  2^16^,1222  +  Pn<f>,2222  =  0,  (2.13) 

which  involves  all  six  independent  elastic  constants.  It  is  assumed  that  these  elastic 
constants  are  such  that  the  associated  strain-energy  density  is  positive-definite  which 
ensures  that  the  partial  differential  equation  (2.13)  is  elliptic  (see  section  2.4  for 
further  discussion).  The  traction-free  boundary  conditions  on  the  surfaces  of  the 
strip  may  be  written  in  terms  of  the  Airy  stress  function  and  then  integrated  to  yield, 
on  using  the  self-equilibration  conditions. 


^  =  0,  ^,2  =  0  at 

X2  —  :tH, 

(2.14) 

4>  =  /(®2),  <f>,i  =  gM 

at  Xi  =  0, 

(2.15) 

0  (uniformly  in  X2) 

as  xi  00, 

(2.16) 

where  f{x2)  and  g{x^  are  prescribed  functions  that  satisfy  certain  continuity  condi¬ 
tions  at  the  corners  of  the  strip.  Since  these  continuity  conditions  will  not  be  used 
here,  their  explicit  forms  are  not  given.  Equations  (2.13)  -  (2.16)  constitute  the 
fundamental  boundary  value  problem  associated  with  the  semi-infinite  anisotropic 
strip  in  a  state  of  plane  stress  or  plane  strain,  which  is  referred  to  henceforth  as  the 
homogeneous  anisotropic  elastic  strip  problem. 
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2.2  Nondimensionalization 


In  the  analysis  of  this  problem  the  governing  partial  differential  equation  is  first  nondi- 
mensionalized,  reducing  the  number  of  independent  elastic  constants.  Returning  to 


equation  (2.13),  we  define  the  dimensionless  coordinates  ?/  as 


(  = 


(2.17) 


and  obtain  the  nondimensional  equation  for  the  Airy  function 


^22 

H* 


2^26 


+ 


(2^12  +  ^66 '\  /^llV 

V  )  W) 


2^16 

H* 


^,nvvri 


=  0. 


(2.18) 


Multiplying  through  both  sides  by  ^  yields 


2^26 


+ 


2^12  +  A 


'66 


L  (Ai/^22)*  J 

2^16 


<l>,Uvv 


^11^22! 


4^,tvvv  d"  ^,vvvv  ~ 


(2.19) 


on  the  domain  ^>0,  — 1<7/<1.  If  one  defines  the  following  nondimensional 
quantities, 


2^12  +  ^66 

-1 

^26 

_1 

016 

.  (^11^22)^  . 

,  ^2  = 

.011022. 

C3  = 

.011022. 

(2.20) 
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then  the  nondimensional  governing  equation  becomes 


2  12 

+  4>,mm  =  0,  (2-21) 

C2  Cl  C3 


on  the  new  region  R  defined  by 

R:  ^  >  0,  -1  <  1/  <  1.  (2.22) 

This  governing  equation  now  involves  the  three  nondimensional  constants  ci,  62  and  €3, 
rather  than  the  six  independent  elastic  constants  ^pq  originally  present.  Transforming 
the  boundary  conditions  (2.14)  -  (2.15)  to  the  nondimensional  variables  results  in 

^  =  0,  =  0  at  Tj  =  ±1,  (2.23) 

4.  =  F(ti),  =  GM  =  0.  (2.24) 

(i>qaP  0  {uniformly  in  X2)  as  ^  00,  (2.25) 

where  a,  ^  £  {^,  j;}  and  F{t]),  G{t])  are  prescribed  functions.  We  note  that  a 
similar  nondimensionalization  scheme  for  an  orthotropic  material  weak  in  shear  was 
carried  out  in  [8].  Similar  nondimensionalizations  have  also  been  used  in  the  buckling 
of  anisotropic  structures  (see  e.g.  [17]  and  the  references  therein). 

2.3  Conservation  Laws 

Before  proceeding  further,  several  inequalities  known  as  “conservation  laws”  are  de¬ 
rived  that  will  be  used  later  to  obtain  the  desired  decay  estimates.  The  basic  method 


CHAPTER  2.  ANISOTROPIC  ELASTIC  STRIP  PROBLEM 


12 


to  obtain  the  inequalities  is  the  same  in  each  case.  The  governing  equation  is  first 
multiplied  by  an  axial  derivative  of  the  Airy  stress  function  and  then  integrated  over 
the  cross  section  of  the  strip.  The  result  is  simplified  using  integration  by  parts  and 
the  known  homogeneous  boundary  conditions  on  the  unloaded  edges  of  the  strip. 

For  the  first  conservation  law,  the  governing  equation  (2.21)  is  multiplied  by 
and  then  integrated  over  where  X*  is  the  cross  section  of  the  nondimensional  strip 
at  ^  =  z,  (see  Figure  2.2),  so  that  we  have 

r  2  12 

/  ~  ~  ~<f>,itivn  +  4>mvv)]  =  0-  (2.26) 

JLz  £2  Cl  ^3 

1 

-1 

Lz 

Figure  2.2:  Cross  section  of  strip 

Performing  integration  by  parts  on  some  of  the  terms,  making  use  of  the  boundary 
conditions  <j}  =  0,  =  0  at  =  ±1,  and  removing  a  ^-derivative  gives  rise  to 

^  MMi  -  ir,  =  0.  (2.27) 
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Upon  integration,  we  obtain  Conservation  Law  1\  i.e., 

[Mm  -  =  Ci,  (2.28) 

where  Ci  is  a  constant.  The  integral  on  the  left  hand  side  of  (2.28)  is  independent  of 
the  axial  coordinate  and  so  is  a  “conserved”  quantity.  Using  similar  procedures, 
the  second  conservation  law  results  from  multiplying  equation  (2.21)  by  and 
integrating  over  the  region  Lz.  Starting  with 

f  2  12 

L  {<f>mi  -  —^mn  +  ~  dr]  =  O  (2.29) 

JLj  62  6i  63 

and  using  integration  by  parts  along  with  the  homogeneous  boundary  conditions,  one 
obtains 

Upon  integration,  Conservation  Law  2  is  obtained;  i.e., 

dv  =  C2,  (2.31) 


where  C2  is  a  constant. 
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The  constants  Ci  and  C2  are  as  yet  unspecified.  However  by  evaluating  the  inte¬ 
grals  as  z  00,  their  signs  may  be  determined  and  the  conservation  laws  (2.28), 
(2.31)  may  be  rewritten.  To  do  this,  we  make  a  minor  additional  assumption  on  the 
behavior  of  a  third  derivative  of  ^  as  ^  >  00,  namely  that 

=  0(1)  {uniformly  in  rj)  as  ^  —y  00.  (2.32) 

It  then  follows  that 

Cl  =0,  Cj  =  -.lim  2  0.  (2-33) 

since  — ♦  0  as  ^  >  00.  Thus,  the  first  and  second  conservation  laws 

imply 


[  +  4>%r,  +  -  L  +  ~^4v\  (2-34) 


and 


—  L  f  (2.35) 

*/jUz  CJ  JLx  t3 


respectively.  These  results  will  be  used  in  later  sections.  It  is  worth  noting  here  that 
the  nondimensional  parameter  C3  (and  so  /3i$)  does  not  appear  explicitly  in  (2.34), 
while  C2  (and  so  ySae)  does  not  appear  explicitly  in  (2.35).  The  details  of  the  derivations 
outlined  above  may  be  found  in  [18]. 
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2.4  Positive-Definiteness  of  the  Strain-Energy 
Density 

It  was  observed  in  section  2.1  that  the  elastic  constants  ^pg  in  the  governing  equation 
(2.13)  are  such  that  the  associated  strain-energy  density  is  positive-definite,  ensuring 
ellipticity  of  the  governing  partial  differential  equation.  This  will  now  be  discussed  in 
more  detail. 

The  total  strain-energy  for  this  boundary  value  problem  is  denoted  by  J^dA 
where  3?  is  the  original  region  occupied  by  the  semi-infinite  strip.  The  strain-energy 
density  W  can  be  expressed  as 


2IV  —  Tap^aP 


(2.36) 


where  upon  using  (2.4),  one  obtains 

21V  =  +  ^22^22  +  ^66^12  d"  2^12'rilT22  +  ^^16TllTi2  -f  2^26'r227l2-  (2.37) 

This  quadratic  form  may  be  written  in  matrix  notation  as  21V  =  t^Bt,  with 


'  Til ' 

^11  ^12  016 

t  = 

T22 

,  B  = 

021  022  026 

.  '^12  . 

.  061  062  066  . 

(2.38) 


where  B  is  the  same  symmetric  matrix  that  appears  in  equation  (2.4)  and  the  super¬ 
script  T  denotes  the  transpose.  The  usual  physical  assumption  that  the  strain-energy 
density  is  positive- definite  ,  i.e.  lV(r)  >  0,  V  r  ^  0  (with  equality  if  and  only  if 
r  =  0)  is  equivalent  to  the  assumption  that  B  is  a  positive-definite  matrix.  Positive¬ 
definiteness  of  a  symmetric  matrix  requires  that  the  diagonal  elements  be  positive 
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and  that  the  determinants  of  the  leading  principal  submatrices  be  positive  as  well. 
Specifically,  this  requires 

/?lls  ^22 »  ^66  >  0>  (2.39) 

Ai^22  “  ^12  >  0,  (2.40) 

^11^22^66  +  2^16^26^12  ~  ^11^26  ■"  ^22^16  ~  ^66^12  ^  (^-^l) 

A  further  requirement  is  that  the  determinant  of  any  submatrix,  including  principal 
minors,  is  positive.  This  yields,  in  addition, 

^22^66  ~  ^26  ^  (2.42) 

Ai^66  —  ^16  >  0-  (2.43) 

For  materials  transversely  isotropic  about  the  xi  axis  (or  for  specially  orthotropic 
materials),  ^le  =  ^26  =  0,  and  conditions  (2.39)  -  (2.43)  reduce  to 


P22,  ^66  >  0> 

(2.44) 

^11^22  “  ^12  >  0, 

(2.45) 

It  will  now  be  shown  that  positive-definiteness  of  the  strain-energy  density  W 
implies  ellipticity  of  the  governing  equation  (2.13).  (Note  that  (2.13)  is  the  original 
version  of  the  governing  equation  rather  than  the  nondimensionalized  version  (2.21). 
When  working  with  the  strain-energy  it  is  more  convenient  to  use  the  original  xi  -  X2 
domain  since  the  strain-energy  is  not  easily  expressed  in  terms  of  the  dimensionless 
parameters  Ci,  €2,  £3  in  the  ^  domain;  see  equation  (2.59)  below.)  The  governing 
partial  differential  equation 
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—  2/326^,1112  +  (2^12  +  ^66)^,1122 

—  2^16^,1222  +  ^11^,2222  =  0  (2.46) 

is  elliptic  provided  that  the  associated  characteristic  equation  has  no  real  roots^  where 
the  characteristic  equation,  denoted  by  P(/i)  =  0,  is 

P(/i)  =  ^ii/i^  +  2/016/1^  +  (2/012  +  +  2/326/1  +  P22  =  0.  (2*47) 

As  noted  in  Lehknitskii  [16]  (see  also  [19]),  choosing  =  [/i^,  1,  nY  and  using  the 
positive-definite  assumptions  on  B  yields 

s^Bs^  =  P{n)  >0  V  /I  real.  (2.48) 

Thus  P{fi)  can  have  no  real  roots  and  (2.46),  i.e.  (2.13),  is  elliptic.  The  estimates  of 
concern  in  this  dissertation  are  relevant  to  elliptic  partial  differential  equations. 

The  positive  definite  conditions  (2.39)  -  (2.43)  will  be  used  in  later  sections  to 
simplify  expressions.  Again  it  is  noted  that  these  conditions  are  in  terms  of  Ppq 
rather  than  ei,  €2  and  €3.  Ideally,  one  would  like  to  obtain  these  conditions  in  terms 
of  the  dimensionless  variables,  but  this  is  not  easily  done  as  the  strain-energy  density 
W  contains  the  constants  /0i2  and  /066  iii  separate  terms  rather  than  combined  as  they 
appear  in  ei  and  the  governing  equation. 

2.5  Energies 

In  this  section,  the  ideas  behind  the  energy  methods  used  in  later  chapters  will  be 
developed.  First,  the  methods  themselves  are  outlined  and  then  several  specific  energy 


functionals  are  derived. 
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2.5.1  Energy  Methods 

The  stress  decay  estimates  presented  in  the  following  chapters  make  use  of  energy 
arguments  that  produce  bounds  on  energy  functionals  and  thereby  lead  to  decay  rate 
estimates  for  the  stresses  in  the  material.  To  motivate  the  choice  of  energy  functionals 
used  in  this  dissertation,  consider  first  an  isotropic  strip  in  a  state  of  plane  stress  or 
plane  strain.  On  using  (2.9),  (2.10)  and  (2.20),  we  find  that 

«.  =  5-  7  =  7  =  »>  (2.49) 

Z  C2  C3 

and  so  (2.21)  reduces  to  the  biharmonic  equation, 

=  0.  (2.50) 

A  natural  energy  functional  associated  with  the  biharmonic  equation  is  the  quadratic 
functional 

■^i(^)  =  /  <l>,ap<f>,ap  dA,  (2.51) 

jRt 

where  the  region  Rz  is  the  shaded  region  illustrated  in  Figure  2.3.  In  (2.51),  the  indices 
a  and  P  denote  either  ^  or  t/  and  the  repeated  index  denotes  summation  over  ^  and  rj. 


Figure  2.3:  Region  Rz 


This  energy  functional  is  related  to  the  total  strain-energy  of  the  material  stored 
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to  the  right  of  ^  =  z  (see  [11])  and  is  referred  to  as  a  “first-order  energy”,  (see  [15]). 
Another  energy  functional  associated  with  the  biharmonic  equation  and  referred  to 
as  a  “second-order  energy”  has  also  been  introduced  (see  [15])  as 

^2(2)  =  /  dA.  (2.52) 

Jr, 

This  energy  can  be  viewed  as  the  first-order  energy  defined  on  The  basic  idea 
behind  the  energy  methods  used  in  later  chapters  of  this  dissertation  is  to  consider  a 
function  F(z)  that  is  a  combination  of  these  energy  functionals.  The  goal  is  then  to 
find  the  largest  positive  value  of  the  constant  k  such  that 

F\z)  -h  2kF{z)  <  0,  (2.53) 

where  the  prime  denotes  differentiation  with  respect  to  z.  This  may  be  accomplished 
through  the  use  of  the  governing  equation,  conservation  laws  and  various  other  in¬ 
equalities.  Equation  (2.53)  implies  the  first-order  differential  inequality 

F'{z)  <  -2kF{z),  (2.54) 

which  leads  to  the  exponential  decay  estimate 

F{z)  <  F{0)e-^^^,  z  >  0.  (2.55) 

It  can  further  be  shown,  depending  on  the  function  F[z),  that  the  energies  themselves 
have  the  same  exponential  behavior,  i.e. 

Eciz)  <  z  >  0,  (2.56) 

where  Ka  (o  =  1,2)  are  constants.  Since  the  energies  are  quadratic  in  the  derivatives 
of  <f>  and  the  stresses  r  are  defined  as  second  derivatives  of  <f>,  one  can  obtain  bounds 
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on  the  stresses  t  such  that  the  exponential  decay  from  the  end  z  =  0  is  of  the  form 

r  ~  Ke-'‘\  (2.57) 

The  value  of  this  estimated  decay  rate  k  gives  a  lower  bound  for  the  decay  rate  of  the 
stresses  in  the  material  occupying  the  domain  in  Figure  2.2. 

2.5.2  Derivation  of  Energy  Functionals 

In  the  previous  section,  the  first-order  and  second-order  energies,  (2.51)  and  (2.52), 
were  described  for  an  isotropic  material.  Their  anadogs  for  the  anisotropic  strip  as 
well  as  additional  energy  functionals  will  now  be  derived. 

There  axe  two  types  of  energy  functionals  that  are  used  in  this  dissertation.  The 
first  type  involves  the  real  physical  strain-energy  of  the  material.  The  second  type 
involves  functionals  that  are  not  directly  related  to  the  physical  strain-energy  but 
are  similar  in  structure.  These  "mathematical  energies”  are  essentially  norms  for 
the  problem.  For  the  anisotropic  strip,  the  strain-energy  density  W  was  expressed 
in  terms  of  the  stresses  t  in  equation  (2.37)  of  section  2.4.  The  total  strain  energy 
E  =  I  WdA  is  written  in  terms  of  the  Airy  stress  function  (2.11)  as 

2E  =  [  2W  dA 

=  +  ^66^fl2  + 

—  2^16^,22^,12  ~  2^6^,11^,12]  dA.  (2.58) 

This  energy  E  is  the  total  physical  energy  of  the  semi-infinite  strip  in  Figure  2.1, 
while  E{z)  is  the  stored  physical  energy  in  the  strip  region  to  the  right  of  X2  =  z. 
Using  the  same  nondimensionalization  scheme  presented  in  section  2.2,  and  letting 
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Ep{z)  now  represent  the  nondimensional  scaled  physical  energy,  one  obtains 

+  'I’m  + 

2  2 

~  ^£v^,riv]  (2.59) 

£2  €3 

where  Rz  is  the  shaded  region  shown  in  Figure  2.3.  The  term  4>,ii<i>,r,Tt  in  the  integral 
(2.59)  may  be  integrated  by  parts  2md  simplified  using  the  boundary  conditions  to 
yield 

/  dy  (2.60) 

JRg  JRx  *'—1 

where  the  last  term  in  (2.60)  is  a  line  integral  contribution.  Thus,  one  obtains 
Hp{z)  =  [<f>%  +  4>\n  +  ~  ~^Ari4>Ai  ~  didr} 

Either  of  the  forms  (2.59)  or  (2.61)  may  be  used  to  express  Ep{z).  In  (2.61),  we 
see  the  explicit  dependence  on  the  parameter  ci,  while  in  (2.59)  this  parameter  is 
not  explicitly  present.  As  was  mentioned  previously,  this  physical  energy  cannot  be 
easily  expressed  solely  in  terms  of  the  parameters  cj  ,£2  smd  £3.  A  similar  observation 
was  made  in  the  context  of  buckling  of  anisotropic  plates  [20].  This  motivates  the 
use  of  mathematical  energy  functionals  that  can  be  expressed  solely  in  terms  of  the 
nondimensional  parameters. 

Next,  several  mathematical  energies  that  are  analogs  to  (2.51),  (2.52)  will  be 
developed.  The  baisic  method  in  developing  these  functionals  is  to  multiply  the  gov¬ 
erning  differential  equation  by  the  Airy  stress  function  itself  or  an  axial  derivative  and 
then  integrate  over  the  region  Rz  in  Figure  2.3.  This  method  is  very  similar  to  the 
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method  used  to  derive  the  conservation  laws  in  section  2.3,  except  that  the  integration 
now  is  over  a  two-dimensional  region  rather  than  a  one-dimensional  region  (line  inte¬ 
gral).  Upon  use  of  the  divergence  theorem  and  simplification  due  to  the  homogeneous 
boundary  conditions,  the  energy  functional  can  also  be  conveniently  expressed  as  a 
line  integral.  The  following  analysis  is  carried  out  in  the  nondimensional  domain. 

To  obtain  the  first-order  energy  functional  for  the  anisotropic  strip,  the  nondi- 
mensionalized  governing  equation  (2.21)  is  multiplied  by  (j>  and  integrated  over  Rg^ 

I  ~  (2.62) 

jRx  C2  ^1  ^3 

Integration  by  parts  and  use  of  the  divergence  theorem  once  yields 

.  2  2  1 

jR,  C2  =3 

r  2  2  1 

=  —  /  (2.63) 

JdR,  C2  ^3  ^3 


where  dRg  denotes  the  boundary  of  Rg,  and  n^,  tifj  are  components  of  the  outward 
unit  normal  to  the  boimdary  in  the  t)  directions.  Repeating  this  process  results  in 


+ 

C2 

ex  ' 

C3 

/  ■ 

+ 

1 

—Mav  - 

1 

JdR, 

C2 

^2 

Cl 

+ 

^3 

^3 

-  Mm  + 

2 

C3 

The  quadratic  functional  on  the  left-hand  side  is  the  anisotropic  analog  of  (2.51).  It 
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will  be  called  a  “first-order”  anisotropic  energy  and  denoted  by  Ei{z),  i.e. 

Ei(z)  =  I  didr}.  (2.65) 

When  Cl  =  |  ,  cj'  =  0  ,  cj ^  =  0  ,  (2.65)  reduces  to  (2.51).  The  boundary 
integral  on  the  right-hand  side  of  (2.64)  may  be  evaluated  and  simplified  using  the 
boundary  conditions  on  the  lateral  sides  of  the  strip  and  at  the  far  end,  so  that  the 
energy  functional  can  be  expressed  in  terms  of  a  line  integral  as 

r  12  2 

Ei{z)  =  -  [-<f><f>m  +  +  —<t>,r,4>Av  +  —HMn  -  dr},  (2.66) 

JLz  Cl  €2  C2 

where  Lz  denotes  the  line  segment  — 1  <  »;  <  1,  i  =  z  (see  Figure  2.2).  Noting 
that 

/  dr)  =  -  f  <i>,n<l>ii  dr)  (2.67) 

JLz  JIjz 

and 

^  dr)  =  [^1  ^dr)  =  0,  (2.68) 

we  get 

Ei{z)  =  -  I  dr).  (2.69) 

JLz  Cl  C2 

It  is  worth  noting  that  the  parameter  €3  does  not  appear  explicitly  in  (2.69).  The 
energy  Ei{z)  given  by  (2.65)  or  (2.69)  is  one  of  two  mathematical  energies  that  will 
be  used  in  later  chapters. 

To  obtain  a  second  anisotropic  energy  functional,  the  governing  equation  is  mul- 
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tiplied  by  and  integrated  over  the  region  R^, 
t  2  12 

I  ~  ~  (2.70) 

jRr  £2  ^3 

Following  the  same  procedure  as  before,  using  integration  by  parts  and  the  divergence 

theorem  repeatedly,  one  obtains 

/  2  12 

J  Rx  ^2  ^3 

JdRx 

2  2  1 

+  - ~  (2*71) 

£2  ^3  ^1 

The  quadratic  functional  on  the  left-hand  side  of  (2.71)  is  the  anisotropic  analog  of 
(2.52).  It  will  be  called  a  “second-order”  anisotropic  energy  and  denoted  by  £^2(^)7 
i.e. 

~  Ir  ~  (2-72) 


The  boundary  integral  on  the  right-hand  side  of  (2.71)  may  be  evaluated  and  simpli¬ 
fied  using  the  boundary  conditions  on  the  lateral  sides  of  the  strip  and  at  the  far  end, 
so  that  the  energy  functional  can  be  expressed  in  terms  of  a  line  integral  as 

E2{z)  =  -  f  (2.73) 

J  L/jf 

We  note  that  none  of  the  nondimensional  material  parameters  ti  (i  =  1,2,3)  ap¬ 
pear  explicitly  in  (2.73).  The  energy  E2{z)  given  by  (2.72)  or  (2.73)  is  the  second 
mathematical  energy  that  will  be  used  in  this  dissertation. 
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It  should  by  noted  that  there  are  many  additional  energy-like  functionals  that  can 
be  obtained  following  the  procedure  used  above.  Depending  on  how  the  integration 
is  performed,  terms  that  are  different  than  those  presented  here  may  arise.  However, 
the  energy  functionals  (2.65),  (2.72)  have  been  found  to  be  the  most  useful  in  the 
present  study.  It  is  also  noted  that  some  of  these  alternate  energy  forms  may  not  be 
positive-definite  and  hence  would  not  be  suitable  norms. 

2.5.3  Positive-Definiteness  of  the  Energy  Functionals 

It  was  stated  previously  that  the  real  physical  strain-energy  of  the  material  must  be 
positive-definite  and  the  conditions  that  this  imposed  on  the  elastic  constants  Ppq 
were  derived.  These  conditions  result  from  physical  properties  and  must  hold  for  all 
materials  regardless  of  which  energy  functionals,  real  or  mathematical,  are  used  in  the 
estimation  techniques.  H  other  energy  functionals  different  from  the  physical  strain 
energy  are  introduced  into  the  problem,  then  this  may  impose  additional  constraints 
on  the  elastic  parameters.  Since  the  mathematical  energies  presented  in  the  previous 
section  are  to  be  used  as  norms,  by  definition  they  must  be  positive.  This  requires 
that  the  integrands  be  positive-definite  quadratic  forms  in  their  arguments.  This  in 
turn  imposes  restrictions  on  Cj,  €2  and  €3. 

Returning  to  (2.65),  one  can  rewrite  the  integrand  of  Ei{z)  as  a  quadratic  form, 
where 


■  1 

0 

1  -1 
\/2C2 

4>i  = 

4>,vri 

and  Bi  = 

0 

1 

- 1 

1 

(2.74) 
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Positive-definiteness  of  Ei  (2)  imposes  the  following  restrictions 

-  >  0,  (2.75) 

Cl 

1  1  1  n 

- 2  -  -2  >  0 

Cl  cl  4 

(We  observe  that  (2.75)  is  implied  by  (2.76).)  Returning  to  (2.72),  the  integrand  of 
E2{z)  can  be  expressed  as  a  quadratic  form,  02®2<I^2?  where 


■  1 

0 

1  1 
N^e2 

<l>2  — 

and  B2  = 

0 

1 

1 

Taia 

72cJ 

1 

>/2e3 

1 

2Ci  J 

(2.77) 


Positive-definiteness  of  E2{z)  again  results  in  the  requirement  (2.76)  since  B2  =  Bi. 

The  restrictions  (2.75)  and  (2.76),  which  result  from  using  the  mathematical  ener¬ 
gies  Ei{z)  and  E2{z),  will  be  assumed  henceforth  and  are  in  addition  to  the  constraints 
(2.39)  -  (2.43)  imposed  by  positive-definiteness  of  the  physical  strain  energy.  All  of 
these  constraints  will  be  useful  in  simplifying  expressions  appearing  in  later  chapters. 


Chapter  3 


THE  ORTHOTROPIC  CASE 


In  the  previous  chapter,  the  elastic  strip  problem  was  formulated  for  anisotropic  ma¬ 
terials.  However,  the  problem  is  simplified  when  one  considers  specially  orthotropic 
and  transversely  isotropic  materials.  These  form  a  special  cleiss  of  materials  for  which 
some  of  the  elastic  constants  vanish.  This  special  case  will  be  examined  first  and  will 
be  referred  to  as  the  orthotropic  case  (although  it  includes  both  specially  orthotropic 
and  transversely  isotropic  materials).  Several  estimation  techniques,  both  analyti¬ 
cal  and  numerical,  will  be  presented.  These  estimates  will  be  compared  with  exact 
numerical  solutions  and  the  results  discussed. 


3.1  Formulation 

The  governing  partial  differential  equation,  conservation  laws,  and  energy  functionals 
of  the  previous  chapter  have  been  derived  for  an  anisotropic  strip.  Consider  now  a 
subclass  of  materials  for  which  the  anisotropic  material  constants  /9i6  and  ^26  vanish. 


27 
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i.e. 

^16  =  1026  =  0.  (3.1) 

(As  pointed  out  in  Chapter  2,  this  simplification  includes  both  specially  orthotropic 
and  transversely  isotropic  materials.)  In  view  of  the  definition  (2.20),  Equation  (3.1) 
implies  that 

i  i  =  0.  (3.2) 

€2  C3 

Returning  to  the  governing  equation  (2.21)  and  using  (3.2),  one  obtains  the  reduced 
equation 

+  ^,vvvv  =  0,  (3.3) 

Cl 

which  contains  only  one  elastic  parameter  €1  instead  of  the  three  parameters  ei,  62, 
C3  present  in  the  anisotropic  case.  (As  we  remarked  in  Chapter  2,  when  ej  =  |  in 
(3.3),  the  familiar  biharmonic  equation  (2.50)  governing  isotropic  materials  results.) 

The  corresponding  reduced  forms  of  the  conservation  laws  may  be  obtained  from 
(2.34)  and  (2.35)  upon  using  (3.2).  This  results  in  the  following  orthotropic  conser¬ 
vation  properties, 


and 


The  strain-energy  density  W  must  be  positive-definite  and  the  restrictions  this  im- 
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poses  on  the  elastic  constants  ySp,  for  the  orthotropic  Ccise  were  presented  in  section 
2.4,  see  equations  (2.44)  and  (2.45). 

The  orthotropic  analogs  of  the  energy  functionals  presented  in  section  2.5  may 
also  be  obtained  by  using  (3.2).  The  nondimensionalized  physical  strain-energy  Ep{z) 
given  by  (2.61)  becomes 

Mil 

where  the  superscript  orth  refers  to  orthotropic  materials.  Again,  it  is  noted  that  this 
physical  energy  cannot  be  written  solely  in  terms  of  the  dimensionless  elastic  param¬ 
eters,  owing  to  the  line  integral  contribution.  The  mathematical  energies  Ei{z)  and 
E2{z),  given  by  (2.65)  and(2.72),  become 

=  L  Wf  + 

=  L. 

respectively.  Using  (2.69)  and  (2.73),  we  can  also  express  these  energies  as  line  inte¬ 
grals  where 

E^*^{z)  =  —  f  [-#,««  +  (3-9) 

JLt  Cl 

Er^{z)  =  -  f  (3'10) 

J  Lig 

The  assumption  (2.75),  i.e.  Ci  >  0,  ensures  that  the  energies  (3.7),  (3.8)  are  positive- 
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definite.  Several  energy  estimates  using  these  energy  functionals  will  be  presented  for 
•  this  orthotropic  elastic  strip  problem.  We  note  that  these  orthotropic  versions  of 

the  governing  equation,  conservation  laws  and  energy  functionals  all  involve  only  one 
elastic  parameter  ti. 


3.2  Preliminary  Inequalities 

Before  proceeding  further,  it  is  convenient  to  record  several  inequalities  which  play 
an  important  role  in  the  arguments  to  follow.  The  first  group  consists  of  the  well 
known  Wirtinger-type  inequalities,  which  have  been  widely  used  in  previous  work  on 
Saint- Venant’s  Principle  (see  e.g.  [1],  [2],  [3]).  These  inequalities  hold  for  sufficiently 
smooth  functions  w(t])  defined  on  the  domain  (—1, 1)  and  are  the  foUowing: 

(i)  li  w{ij)  €  C'^(— 1,1)  and  w{—l)  =  0,  t<;(l)  =  0,  then 

tof,  df)  >  ^  dt].  (3.11) 

(ii)  Ifw(Tf)  €  C'^(— 1,1)  and  u;(— 1)  =  to,,,(— 1)  =  0,  u;(l)  =  t«,»,(l)  =  0,  then 

/I  ri 

(3-12) 

(Hi)  liw{T])  €  C'^(— 1,1)  and  tw(— 1)  =  i«,,,(— 1)  =  0,  ttf(l)  =  «;,,(!)  =  0,  then 


(3.13) 
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where  /Xq  =  4.73,  a  value  slightly  larger  than  For  convenience  this  latter  quantity 
will  be  substituted  in  (Hi)  so  that  one  obtains 

/.'i”™''’'  >  (I)  F 

In  all  of  these  inequalities,  k  is  the  length  of  the  interval  of  integration,  i.e.  h  =  2. 
When  applied  to  the  semi-infinite  strip  problem,  h  represents  the  nondimensionalized 
height  of  the  strip.  These  inequalities  may  be  derived  from  variational  formulations 
and  the  optimal  constants  that  appear  on  the  right-hand  sides  of  (3.11)  -  (3.13)  cor¬ 
respond  to  the  smallest  eigenvalues  in  such  formulations.  Proofs  of  these  inequalities 
may  be  found  in  [21],  [22]  and  [23]. 

Another  inequality  that  will  prove  useful  in  the  subsequent  analysis  is  the  arithmetic- 

geometric  mean  inequality  which  has  the  following  two  forms, 

ta 

—  2a6  <  aa^  -|-  — ,  (a  >  0)  (3.15) 

a 

or 

2ab  <  aa^  -f  — ,  (a  >  0).  (3.16) 

OL 

Both  of  these  forms  may  be  obtained  by  observing  that 

^\/ac  i  (^•^'^) 

3.3  Analytic  Estimate 

The  first  energy  estimate  that  will  be  developed  is  an  analytic  estimate  which  yields 
an  explicit  formula  for  the  estimated  decay  rate  in  terms  of  the  elastic  parameter  Ci. 
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The  estimate  will  be  derived  first  for  the  general  Ccise  of  Ci  >  0  and  then  specialized 
to  an  asymptotic  result  for  small  Ci . 


3.3.1  General  Estimate 


For  this  analytic  estimate  only  one  of  the  energy  functionals,  E^*''{z),  defined  in  (3.7) 
is  used.  Following  an  energy  method  first  developed  in  [11]  for  the  isotropic  case,  let 
the  function  F{z)  be  defined  by 

F{z)  =  Ei{z)  +  2k  Ei{s)ds,  (3.18) 


where  A;  >  0  is  an  unspecified  constant.  Here,  the  superscript  orth  has  been  dropped 
for  convenience  but  it  will  be  assumed  for  the  remainder  of  this  chapter  that  the 
functionals  refer  to  orthotropic  versions.  We  then  have 

F\z)  +  2kF{z)  =  E[{z)  -  2kEi{z)  +  2kEi{z)  +  4P  Ei{s)ds 

=  E[{z)  +  4k^  Ei{s)ds.  (3.19) 

In  order  to  express  the  terms  Ei{z)  and  J  Ei{s)ds  as  line  integrals,  we  first  note  a 
result  derived  in  [11]  that  follows  from  the  definition  of  derivative. 


A. 

dz 


(3.20) 


for  functions  /  continuous  on  Rz.  Using  /  =  “  (3-20) 

yields  an  expression  for  E'-^{z)  as  a  line  integral. 


(3.21) 
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Similarly,  using  /  =  [-<f><f>m  +  (3*20)  implies  in  addition, 

using  (3.9),  that 

Hi{z)  =  —  f  [— (3.22) 
az  JRs  Cl 

On  integration,  and  using  the  divergence  theorem  and  boundary  conditions,  one  ob- 
tains  an  expression  for  J  Ei{s)ds  as  a  line  integral, 

roo  r  1 

/  Ei{s)ds  =  -  [-<l><f>m  + 

Jz  JRt  Cl 

=  -  [Hm  -  d:,.  (3.23) 

Returning  to  (3.19),  we  then  have 
F\z)  +  2kF{z) 


=  - Wm  +  ffi.  + 


,C 


^4%]  dr,.  (3.24) 


Recall  from  section  2.5  that  the  objective  is  to  find  the  largest  positive  value  for  k 
such  that 

F'{z)  +  2kF{z)  <  0.  (3.25) 

Denoting  the  right-hand  side  of  (3.24)  by  J(z)  and  completing  squcires,  one  obtains 

=  -  I  [(<A,«  +  2k^<f>y  + 

JLx 

2k^ 

_  _  4jbV^]  dr]  (3.26) 

which  upon  dropping  the  first  squared  term  yields  the  inequality 

Jiz)  <  -  I  +  4\r,r,  -  dr). 

JLtg  Cl  ti 


(3.27) 


CHAPTER  3.  THE  ORTHOTROPIC  CASE 


34 


CHAPTER  3.  THE  ORTHOTROPIC  CASE 


35 


rate  for  stresses,  we  seek  the  largest  possible  k  in  order  to  yield  the  best  estimate. 
This  reduces  to  the  following  problem  characterizing  the  appropriate  choice  of  k: 


Choose  the  largest  k  such  that 


(3.35) 


subject  to  the  constraint 


‘Ihyjej^ 


(3.36) 


In  these  expressions,  h  =  2  is  the  non-dimensional  strip  width.  Rewriting  (3.35) 


,  l2_?[ _ —  <  0 

*  +  *  2h^e^  h*  - 


(3.37) 


it  is  observed  that  the  largest  value  of  k  for  which  this  is  satisfied  occurs  when  equality 


holds  in  (3.37).  Using  the  quadratic  formula  yields  the  roots  of  this  equation  as 


-  “4^  ('  "  vTTTii) 


(3.38) 


so  that 


5-i)^ 


(3.39) 


where  signs  have  been  chosen  so  that  k  and  P  are  positive  real  quantities.  With  k 
given  by  (3.39),  the  condition  (3.36)  is  satisfied  if  ei  <  When  ei  exceeds  this 


value,  we  choose 


2hy/^' 


f  ^  ^ 

for  Cl  >  . 

4 


(3.40) 


It  is  easily  verified  that  the  choice  of  k  given  by  (3.40)  satisfies  (3.35). 
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From  (3.25),  we  obtain  the  decay  estimate 

F(z)  <  F(0)c-2*^  z  >  0.  (3.41) 

It  will  now  be  shown  that  the  energy  functional  Ei(z)  has  this  same  exponential 
decay  behavior.  Following  the  arguments  of  [11],  it  follows  from  the  definition  (3.18) 
of  F(z)  that 

Ei(z)  <  F(0)e-^’‘\  (3.42) 

It  remains  to  find  an  upper  bound  for  F(0).  Substituting  the  definition  of  F(z) 
directly  into  the  left-hand  side  of  (3.41)  leads  to  the  inequality 

Ei{x)dx^  <  F{0)e-^’‘\  (3.43) 

Integrating  this  inequality  from  2  =  0  to  ir  =  oo  results  in 

Ei{x)dx  <  ^  ^  F;i(0)  -I-  2k  Ei{x)dx  .  (3.44) 

Solving  for  the  integral  of  Ei  gives 

2k  r  Ei{x)dx  <  Ei{0).  (3.45) 

Jo 

Combining  (3.45)  with  (3.18)  yields  the  following  bound  on  F(0), 

F{0)  <  2Ei{0),  (3.46) 

which  when  substituted  back  into  (3.42)  yields  the  desired  result 

Ei{z)  <  2Ei{Q)e-^^%  z  >  0.  (3.47) 

Notice  that  when  z  =  0,  the  inequality  (3.47)  is  not  sharp.  It  has  been  shown  in  [24] 
that  the  multiplicative  factor  of  2  on  the  right-hand  side  of  (3.47)  can  be  eliminated 
(see  also  the  discussion  on  p.  228  of  [1]),  thus  remedying  this  shortcoming. 
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It  is  reasonable  to  expect  from  the  definition  of  Ei(z)  as  quadratic  in  the  second 
derivatives  of  ^  and  hence  quadratic  in  the  stresses  t,  that  the  stresses  themselves 
have  the  exponential  decay  behavior 


kK,  i))|  <  Ke~'^  for  fixed  i), 


(3.48) 


where  r  denotes  a  typical  stress  component.  The  derivation  of  such  pointwise  esti¬ 
mates  is  technically  elaborate  (see  [1]  for  a  discussion)  and  we  shall  not  pursue  this 
here.  The  estimated  decay  rate  is  given  by  (3.39),  (3.40)  as 


k 
k  = 


=  -  0’ 


2hy/t{ 


for  ei  < 
for  Cl  >  — . 


(3.49) 


This  value  of  k  gives  a  lower  bound  on  the  nondimensional  decay  rate  for  stresses 
in  the  material.  Transforming  this  result  back  into  the  original  variables  Xi  and  X2, 
using  the  change  of  variable  presented  in  section  2.2,  implies  that  the  exponential 
decay  factor  in  the  stresses  is  given  by 


(3.50) 


where 


‘‘  =  4^ (V'  +  -  0’  (I;)'  -  X’ 


k*  = 


4-y^  \^22/ 


for  ei  >  ^  (3.51) 


To  obtain  these  expressions,  we  have  set  h  =  2.  We  note  that  ^  «  0.433.  The 
fraction  ^  is  the  estimated  decay  rate  in  the  original  domain  where  H  denotes  the 
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half-width  of  the  original  strip.  The  dependence  of  the  decay  rate  on  the  strip  width 
may  be  more  conveniently  expressed  by  writing  the  decay  rate  as  where  2H  is 
now  the  total  width  of  the  strip.  The  analytic  estimate  (3.51)  gives  a  formula  for  the 
decay  rate  in  terms  of  the  elastic  parameter  cj  and  the  quantity  f  *  which  will  be 
referred  to  as  the  “beta  ratio”. 

For  isotropic  materials  where  Ci  =  (^)«  =  1>  we  obtain  from  (3.51)  that 

k*  _  IT  ^  2.2 

H  ~  2V^ 

Note  that  this  is  precisely  the  same  estimate  that  was  obtained  by  Flavin  [12]  for 
the  biharmonic  equation,  using  a  modification  of  the  earlier  arguments  of  Knowles 
[11].  This  is  also  the  result  obtained  by  Oleinik  and  Yosifian  [13],  [14]  using  alternate 
methods.  The  decay  rate  estimate  given  by  (3.52)  improves  upon  the  estimate  in  [11], 
which  gives  a  decay  rate  of  (^— )^^  2F- 


3.3.2  Asymptotic  Estimate 


The  previous  estimate  was  derived  for  ci  >  0.  If  ci  is  small  enough,  a  Taylor 

expansion  may  be  introduced  yielding  au  asymptotic  version  of  the  analytic  estimate. 
Recalling  the  Taylor  expansion  for  y/1  +  x  as 

Vl  +  ®  ~  ^  2  ~  48^^  +  •  •  •  (l®l  <  1)  (3.53) 

and  substituting  this  expansion  into  the  first  equation  of  (3.51),  with  x  =  Ificj, 
yields 


TT 

4v/^ 


(862  ~  32cj  4-  •  •  •)* 


k* 


(3.54) 
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On  retaining  the  leading  order  term,  we  find  that 


f  jgllV 

V^22/ 


as  £i 


0. 


(3.55) 


Note  that  the  Taylor  expansion  for  k*  is  valid  for  t\  <  It  is  observed  that 
for  small  Cj  this  asymptotic  analytic  estimate  for  k*  approaches  the  general  analytic 
estimate  (3.51)  from  above  due  to  the  alternating  signs  of  the  terms  in  (3.54).  Thus, 
it  is  no  longer  guaranteed  that  (3.55)  is  strictly  a  lower  bound.  In  fact,  if  the  general 
analytic  estimate  is  very  close  to  the  exact  decay  rate,  then  the  asymptotic  estimate 
may  overestimate  the  exact  decay  rate.  Hence,  (3.55)  must  be  viewed  simply  as  am 
approximation  rather  than  a  lower  bound. 

The  values  of  the  two  estimates  (3.51),  (3.55)  for  a  set  of  materials  (in  plane  stress) 
axe  shown  in  Table  3.2.  The  materials  themselves  are  listed  in  Table  3.1  together  with 
their  elastic  moduli.  For  all  the  materials  listed,  the  result  (3.51)i  is  applicable  except 
in  the  isotropic  case,  where  (3.51)2  is  used.  We  will  discuss  these  results  at  the  end 
of  Chapter  3. 


3.4  Conservation  Law  Estimate 

The  next  estimate  that  will  be  developed  is  derived  from  similar  arguments  as  the 
previous  analytic  estimate,  using  in  addition  the  conservation  property  (3.4).  It  will 
be  seen  that  this  yields  a  significantly  improved  result.  Again,  this  estimate  is  first 
derived  for  the  most  general  case  of  ci  >  0  and  then  an  as3anptotic  result  is  derived 
for  smeJl  ci. 
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3.4.1  General  Estimate 

For  this  estimate,  both  of  the  orthotropic  energies  Ei{z)  and  E2{z)  are  used.  Let 
F{z)  be  defined  by 

F{z')  =  [F?2('^)  "I"  ^Fi(^zy\  d"  2k  J  [£(2(5)  d"  (3.56) 

where  m  and  k  are  both  positive  constants  to  be  determined.  The  detailed  procedures 
for  this  estimate  follow  closely  the  arguments  given  in  [3]  for  the  isotropic  strip  where 
ej  =  |.  With  F{z)  defined  in  (3.56),  we  have 

F'{z)  +  2kF{z)  =  E2(z)  d-  mE[{z)  -f  j  [-^2(5)  d-  mEi{s)]  ds.  (3.57) 

too 

We  now  wish  to  express  the  quantities  E2{z)  and  J  E2{s)ds  as  line  integrals,  as  was 
already  done  for  Ei{z)  and  /  Ei{s)ds  in  the  previous  section. 

Following  similar  procedures,  let  /  =  +  J  ”  (S-SO)  ‘o 

obtain  a  line  integral  expression  for  E^^iz), 

E2{z)  =  —  j  d-  —<f>%n  +  (3.58) 

■Jz  Cl 

Letting  /  =  (3-20)  yields 

E2{z)  =  ^  ^  (3-59) 

which,  when  combined  with  the  divergence  theorem  and  boundary  conditions,  yields 
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fOO 

an  expression  for  j  E2{s)ds  as  a  line  integral, 

I  E2{s)ds  =  “  /  ~ 

J  Z  J  it* 

=  L  f  1 


(3.60) 


Substituting  (3.21),  (3.23),  (3.58)  eind(3.60)  into  (3.57)  results  in 

f'(z)  +  2kF{z)  =  -  /  +  <l>\„  +  (to  - 

J  Lx 

+  (m  -  2P)<^f,,  + 

=  -J{z).  (3.61) 


We  now  seek  positive  values  for  k  and  m  such  that  the  integral  J(2)  is  nonnegative. 
The  first  step  in  simplifying  J{z)  is  to  make  use  of  the  conservation  property  (3.4) 
which  may  be  rewritten,  upon  using  the  arithmetic-geometric  mean  inequality  (3.15), 
as 

/  [<f>%v  -  ~  -  ~Jr 

JLx  Cj  j  Lz 

<  rf,  (3.62) 

for  a  >  0.  Thus  we  have  from  (3.62)  that 

L  -  L 

Using  (3.63)  in  (3.61),  we  obtain 

J{z)  >  I  [{m-2k^  +  a)<l>%^  +  (m  ~  2P  - 

JLx  ^1 

-  (4fc*m  +  +  4k^m(l><f>,(i 

^1 


(3.64) 


CHAPTER  3.  THE  ORTHOTROPIC  CASE 


42 


On  using  again  the  arithmetic-geometric  mean  inequality  (3.15),  we  obtain 

/  dri  >  -7  /  I  <l>%  dr),  7  >  0. 

J  Lx  JLx  I  *'Lx 


(3.65) 


On  employing  the  Wirtinger  inequality  (3.14)  with  the  choice  w  =  <f>,  we  obtain 
from  (3.65) 

ir,  >  <i,  -  i  ■!>%  d^,  (3.66) 


which  may  be  used  in  (3.64)  to  yield 


2k^m, 


Jiz)  >  -I-  +  {m-2k'^-a-^^)(l>\^ 

/  r.,-1  2Pm'yh*..o  .  (m  —  a),2 

+  (m-2e  +  a- 

-  (4l“m  +  di). 


(3.67) 


Next  several  Wirtinger  inequalities  will  be  used,  the  validity  of  which  make  use  of 
the  boundary  conditions  (2.23).  The  choices  w  =  (f>  and  w  =  in  the  inequality 
(3.12)  yield 


and 


(3.68) 


(3.69) 


respectively.  Similarly,  the  choices  of  lu  = 


in  (3.11)  and  w  =  <f>^^  in  (3.14)  )deld 


(3.70) 
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and 


(3.71) 


respectively.  Using  (3.68)  -  (3.71)  in  (3.67),  provided  that  (m  —  a)  <  0,  we  obtain 


2 

f  _,2  2Pm-i  ,2 

+  ihiT  +  m  -  2k  -  a - —}ht 


f  «t2  2k^m'yh*  Pmh^'t  -i  /« ^o\ 

+  {m  -  2A:2  +  a - - 

It  then  follows  that  the  integral  J{z)  will  be  nonnegative  if  we  can  choose  positive 
constants  m,  a,  7  and  k,  such  that  the  coefficients  of  and  in  (3.72)  are 

nonnegative.  Following  the  notation  introduced  in  [3],  we  let 


TT' 


m  = 


,2  ^^1^2  a  c 

‘  ’  “  =  F^’  ^  =  F^’ 


(3.73) 


where  h  =  2.  The  problem  is  then  to  find  positive  constants  M,  /3,  6  and  K,  such 
that 


1  + 

-  ■+ 
ei 


4ei 


>  M, 

(3.74) 

(IF  " 

(1)^  - 

(3.75) 

2K'^M 

-  /3  - 

S 

(3.76) 

K^M 

2K^M6 

mmm. 

2ci 

■  (1)^ 

(3.77) 

M  -  2K^  +  ^  - 

We  desire  the  choice  of  constants,  satisfying  (3.74)  -  (3.77),  to  maximize  the  value  of 
K  and  hence  maximize  k.  In  [3]  it  was  shown  that  choosing 


^  =  M 


(3.78) 
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resulted  in  the  maocimum  value  of  k  for  the  case  where  ei  =  i.e.  for  isotropic 
materials.  Similar  arguments  can  be  used  to  show  that  (3.78)  also  provides  the 
maximutn  value  of  k  in  the  present  case.  On  choosing  (3.78),  the  remaining  conditions 
(3.75)  -  (3.77)  become 


VI 

Kir  - 

AM 

(3.79) 

VI 

1  r. 

(3.80) 

2ci  + 

2M 

(3.81) 

_ 

-  L 

M 

2M6]  ■ 

+ 


261  ■  (1)^ 


It  is  readily  verified  that  the  right-hand  sides  of  (3.79), (3. 80)  and  (3.81)  are  strictly 
monotone  functions  of  the  variables  M  and  6.  Thus  the  maximum  choice  for  K 
satisfying  (3.79)  -  (3.81)  results  from  equating  these  three  expressions.  Equating 
(3.79)  and  (3.80)  results  in  the  following  expression  for  5, 

2€xM[{\Y  -  M2] 


6  = 


{2ei[(|)4-M2]-4M}- 


(3.82) 


Equating  (3.79)  and  (3.81),  and  using  (3.82),  results  in  a  sixth  order  polynomial 
equation  for  M, 

+  20e.M*  -  (|)V  -  4(|)‘e,]  +  4M  -  2  (|)%.] 

-  [(j)‘  -  tA  =  0.  (3.83) 

The  roots  of  (3.83)  may  then  be  used  to  obtain  K  by  taking  equality  in  (3.79),  where 
it  can  be  observed  that  the  smallest  value  of  M  yields  the  largest  value  for  K.  The 
sixth  order  polynomial  (3.83)  was  solved  using  MACSYMA  for  the  set  of  materials 
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given  in  Table  3.1.  For  each  specific  material,  ci  was  determined  and  substituted  into 
(3.83),  which  was  then  solved  numerically.  For  all  of  the  materials,  the  six  roots  Mi, 
i  €  {1,2, 3, 4, 5, 6}  were  such  that  the  smallest  positive  root  resulted  in  a  negative 
choice  for  S,  and  hence  was  inadmissible.  The  second  positive  root  in  each  case  gave 
the  desired  maximum  value  for  K  satisfying  (3.79)  -  (3.81)  and  hence  from  (3.73)  the 
maximum  value  for  k. 

We  see  from  (3.61)  that  this  value  for  k  yields  an  exponential  decay  for  F{z)  of 
the  form 

F{z)  <  F(0)e-2*^  z  >  0.  (3.84) 

From  the  definition  of  F(z)  in  (3.56),  we  find  that 

Eiiz)  +  mEi{z)  <  F(0)e-=**%  ^  >  0.  (3.85) 

Adapting  the  arguments  given  in  (3.42)  -  (3.46)  and  modelled  after  [11],  it  can  be 
shown  that 

F{0)  <  2[F;2(0)  +  mF;i(0)].  (3.86) 

It  then  follows  from  (3.85)  and  (3.86)  that 

Eiiz)  +  mEi{z)  <  2[JS2(0)  +  mEi{0)]e-^'’%  z  >  0,  (3.87) 

and  hence,  since  Ei{z)  and  E2{z)  are  nonnegative, 

Ei{z)  <  2  [Ei(0)  +  ^^2(0)]  e-^*^  z  >  0,  (3.88) 

—2kz 


E2{z)  <  2[F?2(0)  +  mF?i(0)]e 


5 


z  >  0, 


(3.89) 
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showing  that  the  energies  Ei(z)  and  E2(z)  have  the  same  exponential  decay  as  F{z). 
Again,  it  can  be  shown  using  arguments  in  [24]  that  the  multiplicative  factors  of  2  on 
the  right-hand  sides  of  (3.88)  and  (3.89)  can  be  eliminated. 

The  derivation  of  pointwise  estimates  for  the  stresses  r  of  the  form  (3.48)  is  again 
technically  elaborate  (see  e.g.  [1],  [2],  [3]).  Assuming  a  result  of  the  form  (3.48),  and 
transforming  this  result  back  into  the  original  domain  leads  to  an  estimated  stress 
decay  rate  of  where  k*  is  given  by 

r  =  .(|l)\  (3.90) 

Although  this  conservation  law  estimate  does  not  give  an  explicit  formula  for  the 
decay  rate  in  terms  of  the  elastic  constants,  the  numerical  results  show  a  considerable 
improvement  over  the  analytic  estimate  (3.51).  These  results  are  shown  in  Table  3.2. 
The  results  for  the  isotropic  case,  which  coincide  with  [3],  are  also  shown  in  Table 
3.2. 

3.4.2  Asymptotic  Estimate 

The  technique  based  on  the  conservation  property  (3.4)  can  also  be  used  to  obtain 
an  asymptotic  estimate  for  small  ci.  Using  perturbation  methods,  we  look  for  roots 
M  of  (3.83)  of  the  form 

M  =  mo  +  mi€i  -|-  m2€i  + .  .  .  (3.91) 

where  M  has  been  expanded  in  powers  of  the  small  parameter  ci  and  mo,  mi  and 
m2  are  constants  to  be  determined.  Keeping  only  the  first  three  terms,  we  substitute 
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(3.91)  into  (3.83)  aaid  collect  like  powers  of  Ci.  This  results  in  an  equation  of  the  form 

Ao  +  Aici  +  A2€l  +  ...  +  0(ef+^)  =  0.  (3.92) 

Using  the  Fundamental  Theorem  of  Perturbation  Theory  (see  [25]),  we  set  each  of 
the  coefficients  A,  in  (3.92)  equal  to  zero,  which  in  turn  gives  values  for  mo,  mi  and 


m2.  This  process  results  in  four  roots  of  the  form  (3.91), 

A/,  =  (|)“  -  8£,  +  +  .  .  .  (3.93) 

M,  =  -  (I)'  -  8£.  -  +  .  .  .  (3.94) 

M,  =  1(1)%.  +  ...  (3.95) 

M4  =  -4ci  +  .  .  .  ,  (3.96) 


while  the  remaining  two  roots  of  (3.83)  take  on  a  different  expansion  that  is  propor¬ 
tional  to  It  is  observed  that  Mi  and  M3  give  two  positive  roots  while  M2  and 
M4  give  two  negative  roots.  Furthermore,  it  is  seen  that  Mi  given  by  (3.93)  yields 
the  second  positive  root,  which  vras  the  root  of  (3.83)  obtained  using  MACSYMA 
that  yielded  the  maximum  value  for  K.  Taking  the  first  two  terms  in  (3.93)  as  an 
approximation  for  M,  i.e. 

M  ~  7  —  8ei,  as  ei  -4  0,  (3.97) 

4 

cind  substituting  into  (3.79)  with  equality  holding,  we  find  that 

[(1)^  -  (!  -  8ei)^]^ 

2(|-8ci)i 
(36€i  -  64ef)» 

(9  -  32ei)3 

-  f‘i)'  (1  -  f'l)"'’ 
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Noting  the  Taylor  expansions 

(1  +  =  1  +  I  -  y  +  .  . .  (1*1  <  1)  (3.99) 

(1  +  x)-i  =  1  -  I  +  ^  +  . .  .  (1*1  <  1)  (3.100) 

we  obtain  from  (3.98), 

8  16 

K  ~  2Vir(l  -  gCi  +  .  .  .  )(1  +  yei  +  .  .  . ) 

~  2y/^{\  +  -€l  +  .  .  .  ) 

2'y^i7j  9,  (3.101) 


where  the  leading  order  term  has  been  retained  in  the  last  step. 

Returning  to  (3.73)  and  (3.90)  we  have  the  following  asymptotic  conservation  law 
estimate  for  k*, 


k* 


as  €i  — *■  0 


(3.102) 


which  improves  upon  the  asymptotic  analytic  estimate  (3.55)  by  a  factor  of  v^- 
Note  that  this  estimate  is  valid  for  £i  <  ^  ^  -28.  It  has  also  been  shown  in 
[6],  [7]  that  (3.102)  follows  from  an  asymptotic  analysis  of  the  exact  decay  rate,  so 
that  the  conservation  law  argument  yields  an  estimate  for  the  decay  rate  which  is 
asymptotically  exact.  This  asymptotic  conservation  law  estimate  is  shown  for  a  set  of 
orthotropic  materials  in  Table  3.2.  We  again  observe  that  (3.102)  is  not  guaranteed 
to  be  a  lower  bound  for  the  exact  decay  rate  but  rather  an  approximation,  valid  for 
small  Cl.  Further  discussion  of  the  result  (3.102)  will  be  given  later  in  this  section. 
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3.5  Exact  Decay  Rates 

The  exact  solution  to  (3.3)  subject  to  the  boundary  conditions  (2.23)  -  (2.25)  may 
be  obtained  using  a  separation  of  variables  technique.  First,  we  seek  solutions  of  the 
form 

<f)  =  7  =  constant.  (3.103) 

This  form  is  chosen  since  we  seek  stresses  that  decay  exponentially  in  the  axial  direc¬ 
tion.  Substituting  this  form  for  <f>  directly  into  the  governing  equation  and  boundary 
conditions  gives 

F""  +  iLp"  +  =  0,  on  (-1,1) 

Cl 

F{-1)  =  F(l)  =  0,  (3.104) 

F'i-1)  =  F'(1)  =  0. 

This  is  a  fourth-order  eigenvalue  problem  for  an  ordinary  differential  equation.  Solu¬ 
tions  to  (3.104)  are  sought  in  the  form 

F  =  Ae”*’’,  m  =  constant,  (3.105) 

which  when  substituted  back  into  the  boundary  conditions  (3.104)  yields  one  of  three 
separate  characteristic  equations  (or  eigenconditions)  for  7  ,depending  on  the  range 
of  t\  (see  [6],  [19]).  The  eigenvalues  7  are  complex  in  general.  Numerical  solutions  to 
these  eigenconditions  have  been  computed  in  [6]  and  [19].  The  decay  rate  is  obtained 
from  the  eigenvalue  with  the  smallest  positive  real  part.  Once  obtained,  this  result 
must  be  transformed  back  to  the  original  variables  xi  and  X2  to  arrive  at  the  exact 
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decay  rate  for  the  material.  Exact  decay  rates  for  a  set  of  orthotropic  materials  (in 
plane  stress)  are  shown  in  Table  3.2  using  the  methods  of  [19]  which  will  be  discussed 
more  fully  in  the  next  chapter.  In  this  table,  note  that  for  all  the  materials  other 
than  isotropic,  the  exact  decay  rate  is  just  the  purely  real  first  eigenvalue.  In  fact  it 
can  be  shown  from  [19]  that  the  transition  point  where  the  eigenvalue  changes  from  a 
purely  real  value  to  a  complex  value  occurs  when  Fbr  ^  condition 

satisfied  by  all  of  the  orthotropic  materials  in  the  table),  the  first  eigenvalue  is  purely 
real. 

In  the  case  of  isotropic  materials  where  cj  =  |,  the  eigenvalue  problem  (3.104) 
reduces  to 

F""  +  +  7^F  =  0,  on  (-1,1) 

F(-l)  =  F(l)  =  0,  (3.106) 

F'(-l)  =  ^'(1)  =  0- 

This  is  a  well-known  eigenvalue  problem  whose  solutions  F  are  known  as  the  Fadle- 
Papkovich  eigenfunctions  (see  e.g.  p.  231  of  [1]  for  a  discussion).  Seeking  solutions  of 
the  form  (3.105)  and  using  the  boundary  conditions  (3.106)  leads  to  an  eigencondition 
for  7, 

sin  (27)  =  ±27.  (3.107) 

The  smallest  real  part  of  an  eigenvalue  7  arises  from  the  taking  the  minus  sign  in 
(3.107),  (corresponding  to  symmetric  deformations),  and  is  given  by 


Fe  7  «  2.106, 


(3.108) 
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so  that  the  stresses  decay  as 

(3.109) 

Transforming  this  result  back  into  the  original  variables,  we  see  that  the  stresses  decay 
as 

(3.110) 

and  so  we  recover  the  well-known  result  (see  e.g.  [1],  [26])  that  the  exact  decay  rate 
for  stresses  in  the  isotropic  material  is 

3.6  Discussion  of  Results 

We  now  discuss  the  stress  decay  estimates  for  orthotropic  strips  obtained  in  this 
chapter  and  compare  them  with  the  corresponding  exact  decay  rates  that  were  com¬ 
puted  numerically.  The  materials  used  for  comparison  are  contemporary  engineering 
materials,  and  are  listed  in  Table  3.1  along  with  their  elastic  moduli  and  Poisson’s 
ratios.  All  of  the  materials  are  specially  orthotropic  with  the  fibers  parallel  to  the 
Xi  direction.  Decay  estimates  were  evaluated  for  each  of  the  materials  in  Table  3.1 
and  the  results  axe  shown  in  Table  3.2  along  with  the  exact  decay  rates.  All  of  the 
results  were  computed  for  plane  stress  (and  thus,  from  (2.7),  the  values  of  i/tt  were 
not  needed).  The  isotropic  case  is  included  in  the  last  row  of  the  table  for  comparison. 

The  ansdytic  estimate  (3.51),  asymptotic  analytic  estimate  (3.55)  and  asymptotic 
conservation  law  estimate  (3.102)  all  yield  an  explicit  formula  for  the  decay  rate  in 
terms  of  the  parameter  cj  and  the  beta  ratio  (^)‘ •  (Recall  that  the  general  conser¬ 
vation  law  estimate  did  not  yield  an  explicit  formula  and  thus  was  purely  a  numerical 
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estimate.)  It  is  now  useful  to  express  these  decay  estimates  in  terms  of  the  engineer¬ 
ing  constants  of  the  material.  Using  the  definition  (2.20)  of  Ci  and  recalling  (2.7),  we 
obtain  that  for  plane  stress, 

Here  E,  i/  and  G  are  the  Young’s  Modulus,  Poisson  ratio  and  shear  modulus  respec¬ 
tively.  Using  (3.111),  the  analytic  estimate  (3.51)  becomes 


k*  =  ^ 


TT  /  Et  2Etvlt 


4  V  ^LT  El 


/ 1  + 


El 


Et 

if 


(  1  "7 

\El  —  j 

V  J?T  KEl  —  2vijtGlt'  ~  4  ’ 


r  = 


/  Et 

2Eti^lt 

JGlt 

El  '• 

if 


El  /  Glt  ^  ^ 

Et  \El  —  2vltGlt)  4 

(3.112) 


Similarly,  the  asymptotic  analytic  estimate  (3.55)  becomes 


JL  / 

y/2  V  El  —  2vlt 


2ultGlt  (^^L-21/LTGLr) 


while  the  asymptotic  conservation  law  estimate  (3.102)  becomes 


,,  /  Glt 

e  ^  ’r\/-= - - — 

V  El  —  2ult 
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To  discuss  these  results,  we  refer  to  Table  3.2  where  it  is  observed  from  examining 
the  first  column  that  the  materials  axe  ordered  according  to  increasing  values  of  ci.  It 
is  also  observed  that  all  of  the  decay  estimates  presented  in  the  table  have  the  same 
relative  ordering  as  the  exact  decay  rates  given  in  the  far  right  column.  The  second 
column  shows  the  analytic  estimate  (3.51)  (or  equivalently  (3.112))  for  each  of  the 
materials.  For  the  isotropic  case,  the  analytic  decay  estimate  is  roughly  half  of  the 
exact  decay  rate.  However,  as  one  moves  upward  in  the  table  (i.e.  as  Ci  decreases)  the 
anedytic  estimate  improves  in  accuracy  so  that  for  boron  epoxy  which  has  the  smallest 
value  of  Cl ,  the  analytic  decay  estimate  is  about  70%  of  the  exact  decay  rate.  As  was 
mentioned  previously,  (3.51)i  was  used  for  all  materials  except  for  the  isotropic  case 
where  (3.51)2  Weis  used. 

The  asymptotic  analytic  estimate  (3.55)  (or  equivalently  (3.113))  is  presented  in 
the  third  column  of  Table  3.2.  For  the  materials  for  which  this  estimate  is  applicable 
(i.e.  Cl  <  |),  we  see  a  slight  improvement  over  the  general  analytic  estimate.  Again, 
we  see  that  as  Ci  decreases  the  agreement  between  the  asymptotic  analytic  estimate 
and  the  exact  decay  rate  improves  as  expected. 

The  fourth  column  in  Table  3.2  shows  the  general  conservation  law  estimate  that 
arises  from  solving  the  sixth  order  polynomial  given  by  (3.83).  While  this  estimate 
does  not  provide  an  explicit  formula  for  the  decay  rate,  obtaining  it  is  not  compu¬ 
tationally  intensive.  The  results  show  much  better  agreement  with  the  exact  decay 
rates  than  obtained  by  either  of  the  analytic  estimates.  The  improvement  for  the 
isotropic  case  is  only  slight,  but  as  Ci  decreases  the  improved  agreement  becomes  much 
better.  For  materials  with  small  values  of  Cj  the  conservation  law  estimate  yields  a 
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fairly  accurate  approximation  to  the  exact  decay  rate. 

Finally,  the  fifth  column  shows  the  as5rmptotic  conservation  law  estimate  (3.102) 
(or  equivalently  (3.114)).  For  the  range  of  Ci  where  this  decay  estimate  is  valid  i.e. 
Cl  <  0.28,  we  see  a  further  improvement  over  the  general  conservation  law  estimate. 
Again,  as  ci  decreases  this  asymptotic  estimate  more  closely  approximates  the  exact 
decay  rate  so  that  for  small  Ci  there  is  very  little  error.  This  is  to  be  expected 
since,  as  mentioned  previously,  (3.102)  also  results  from  an  asymptotic  analysis  of  the 
eigencondition  characterizing  the  exact  decay  rate. 

The  results  in  Table  3.2  indicate  that  the  orthotropic  decay  estimates  obtained 
in  this  chapter  yield  the  best  agreement  with  corresponding  exact  results  for  the 
materials  with  small  values  of  ci .  This  may  be  related  to  the  fact  that  for  small  values 
of  Cl  the  smallest  eigenvalue  of  equation  (3.104)  is  purely  real,  and  the  methods  used 
in  this  chapter  make  use  of  real  analysis  arguments  to  estimate  this  quantity.  It  is  also 
observed  that  the  conservation  law  approach  yields  the  best  estimate,  asymptotically 
approaching  the  exact  decay  rate.  This  suggests  that  the  conservation  laws  derived  in 
Chapter  2  play  an  important  role  in  estimating  the  exact  decay  rate  since  they  provide 
results  that  are  in  good  agreement  with  corresponding  exact  results  and  are  not 
computationally  intensive.  For  both  the  analytic  and  conservation  law  estimates,  their 
asymptotic  versions  yield  improvements  suggesting  again  that  the  energy  methods  of 
this  chapter  work  best  for  materials  with  small  values  of  Ci. 

While  it  appears  from  Table  3.2  that  the  decay  estimates  are  ordered  in  terms  of 
increasing  ci  (with  the  boron  epoxy  as  an  exception),  this  is  somewhat  misleading. 
From  the  formulae  (3.51),  (3.55)  and  (3.102),  we  see  that  the  decay  estimates  are 
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dependent  on  both  Ci  and  the  beta  ratio  and  are  thus  ordered  by  a  combination 

of  these  quantities.  An  interesting  simplification  results  when  one  considers  specially 
orthotropic  (or  transversely  isotropic)  materials  with  a  high  degree  of  orthotropy  in 
the  axial  direction.  For  such  strongly  orthotropic  materials,  we  assmne  (see  [6],  [7]) 


that 
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Glt 

El 


«  1, 


Glt 

Et 


<  1. 


Using  (3.115)  in  (3.111)  we  obtain 

^  Glt  { ^11  \  * 

■s/ElEj^  \^22/  \El) 


(3.115) 


(3.116) 


It  follows  from  (3.115)  and  (3.116)  that  ei  «  1  and  that  the  decay  estimates 

for  small  ei  are  valid  for  strongly  orthotropic  materials.  Upon  using  (3.116),  the 
asymptotic  analytic  estimate  (3.55)  and  asymptotic  conservation  law  estimate  (3.102) 
reduce  to 
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JL. 

V2\El)  ’ 
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(3.117) 


and 


as 


Glt 

El 


0, 


(3.118) 


respectively.  Both  of  these  decay  estimates  are  of  the  form 


(3.119) 


a  result  which  Wcis  first  established  in  [5]  using  the  energy  decay  arguments  of  [4]. 
The  result  (3.118)  was  also  obtained  in  [6],  [7]  from  an  asymptotic  analysis  of  the 
exact  decay  rate. 
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The  material  constant  ratios  given  in  (3.115),  the  asymptotic  conservation  law  es¬ 
timate  for  strongly  orthotropic  materials  (3.118),  and  the  exact  decay  rates  (repeated 
from  Table  3.2  for  convenience)  are  tabulated  in  Table  3.3.  In  this  table,  the  mate¬ 
rials  are  now  ordered  by  increasing  values  of  thus  switching  the  order  of  boron 
epoxy  and  ultra-high  modulus  graphite  epoxy  from  Table  3.2.  The  estimate  (3.118) 
and  the  exact  decay  rates  are  seen  to  increase  with  All  but  the  last  two  mate¬ 
rials  in  the  table  may  be  considered  strongly  orthotropic.  The  result  (3.118)  provides 
an  extremely  accurate  estimate  for  the  exact  decay  rate  for  the  strongly  orthotropic 
materials.  (See  e.g.  [2],  [7],  [9],  [10]  for  a  discussion  of  the  utility  of  (3.118)  in  the 
mechanics  of  laminated  composite  structures.) 

The  results  of  this  chapter  have  been  obtained  for  the  case  of  plane  stress.  Nu¬ 
merical  results  for  the  case  of  plane  strain  have  also  been  calculated.  The  results  are 
not  presented  here  as  they  differ  very  little  from  the  plane  stress  case.  Recall  that  the 
only  difference  between  these  two  formulations  is  in  the  values  of  the  elastic  constants 
^pq.  For  the  case  of  strongly  orthotropic  materials,  we  still  obtain  (3.118)  in  plane 
strain. 

The  decay  rate  results  presented  here  may  also  be  expressed  in  terms  of  decay 
lengths.  If  we  denote  the  characteristic  decay  length  d  as  the  length  over  which  the 
solution  <f>  and  hence  the  stresses  r  decay  to  1%  of  their  original  value,  i.e. 

|r(d,a:2)|  =  Ke~^*^  =  for  fixed  X2,  (3.120) 


then  we  obtain  that 
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where  2H  is  one  strip  width  of  the  material.  Thus  the  decay  length  is  inversely 
proportional  to  the  decay  rate  k*,  and  underestimates  of  k*  provide  conservative 
overestimates  of  d.  Table  3.4  shows  the  exact  characteristic  decay  lengths  for  each 
of  the  materials  previously  mentioned.  (The  exaet  decay  rates  are  repeated  from 
Table  3.2  for  convenience.)  From  the  last  row  of  Table  3.4  we  recover  the  well  known 
result  that  the  isotropic  material  has  a  characteristic  decay  length  of  roughly  one 
strip  width,  meaning  that  approximately  one  strip  width  away  from  the  edge  the  end 
effects  are  negligible.  This  is  the  general  rule  often  invoked  when  using  Saint- Venant’s 
principle  for  isotropic  materials  (see  e.g.  [1]).  However,  for  orthotropic  materials  the 
decay  length  may  be  much  larger.  As  seen  in  Table  3.4,  the  decay  lengths  vary 
from  roughly  one  strip  width  to  six  strip  widths  in  the  case  of  ultra-high  modulus 
graphite  epoxy.  For  this  situation,  end  effects  may  persist  much  farther  into  the 
strip  than  for  the  corresponding  isotropic  case.  Thus  for  orthotropic  materials  Saint- 
Venant’s  principle  cannot  be  routinely  applied  as  for  isotropic  materials.  This  fact 
poses  a  serious  deficiency  in  the  use  of  Saint- Venant’s  principle  to  obtain  approximate 
solutions  to  anisotropic  boundary- value  problems.  This  emphasizes  the  importance, 
particularly  in  the  case  of  strongly  orthotropic  materials,  of  obtaining  accurate  decay 
rate  estimates  such  as  (3.118).  Here,  very  high  Mpect  ratios  (ratio  of  length  to  width) 
are  necessary  before  end  effects  can  be  neglected. 
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Material 

El  (psi) 

Et  (psi) 

^LT 

I/XT 

Glt  (psi) 

Boron  Epoxy 

30.0x10® 

3.0x10® 

.21 

.35 

0.7x10® 

Ultra-high  Modulus 
Graphite  Epoxy 

45.0x10® 

0.9x10® 

.26 

NA 

0.6x10^ 

Kevlar  Epoxy 

11.0x10® 

0.8x10® 

.34 

NA 

0.3x10® 

High  Strength 
Graphite  Epoxy  1 

20.0x10^ 

1.0x10® 

.25 

.25 

0.6x10® 

High  Strength 
Graphite  Epoxy  2 

18.5x10® 

1.6x10® 

.35 

.20 

0.832x10® 

S-glass  Epoxy 

7.5x10® 

1.7x10® 

.25 

NA 

0.8x10® 

Boron  Aluminum 

33.0x10^ 

21.0x10® 

.23 

.30 

7.0x10® 

Table  3.1:  Orthotropic  materials  and  elastic  moduli 
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Material 

(«i) 

Analytic 

Estimate 

(3.51) 

Analytic 
Estimate 
for  small 
(3.55) 

Cons.  Law 
Estimate 

Cons.  Law 
Estimate 
for  small 
(3.102) 

Exact 
Decay 
Rate  k* 

BE 

(0.0745) 

0.337 

0.341 

0.445 

0.482 

0.487 

UM 

(0.0949) 

0.253 

0.257 

0.328 

0.364 

0.369 

KE 

(0.1030) 

0.363 

0.370 

0.467 

0.524 

0.534 

HSl 

(0.1362) 

0.375 

0.388 

0.468 

0.548 

0.567 

HS2 

(0.1578) 

0.458 

0.478 

0.564 

0.677 

0.708 

SE 

(0.2366) 

0.684 

0.745 

0.794 

1.054 

1.195 

BA 

(0.2946) 

0.955 

NA 

1.069 

NA 

1.942 

Isotropic 

(0.500) 

1.111 

NA 

1.220 

NA 

2.106 

Table  3.2:  Plane  stress  results  for  orthotropic  materials  :  r  ~  e  as  xi  oo. 


Key  for  Materials 


BE  = 

boron  epoxy 

UM  = 

ultra-high  modulus  graphite  epoxy 

KE  = 

kevlax  epoxy 

HSl  = 

high  strength  graphite  epoxy  1 

HS2  = 

high  strength  graphite  epoxy  2 

SE  = 

s-glass  epoxy 

BA  = 

boron  aluminum 
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Material 

(^0 

Glt/El 

Ex/El 

Glt/Et 

Cons.  Law 
Estimate 
(3.118) 

Exact 
Decay 
Rate  k* 

UM 

(0.0949) 

0.013 

0.02 

0.66 

0.358 

0.369 

BE 

(0.0745) 

0.023 

0.10 

0.23 

0.476 

0.487 

KE 

(0.1030) 

0.027 

0.072 

0.37 

0.516 

0.534 

HSl 

(0.1362) 

0.03 

0.05 

0.6 

0.544 

0.567 

0.045 

0.086 

0.52 

0.666 

0.708 

SE 

(0.2366) 

0.106 

0.226 

0.47 

NA 

1.195 

BA 

(0.2946) 

0.212 

0.636 

0.33 

NA 

1.942 

Table  3.3:  Strongly  orthotropic  results:  r  ~  e  h  as  xi  —y  oo. 
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Material 

(^0 

Exact 
Decay 
Rate  k* 

Exact 

Characteristic 
Decay  Length  d 

UM 

(0.0949) 

0.369 

6.24  X  2H 

BE 

(0.0745) 

0.487 

4.73  X  2H 

KE 

(0.1030) 

0.534 

4.31  X  2H 

HSl 

(0.1362) 

0.567 

4.06  X  2H 

HS2 

(0.1578) 

0.708 

3.25  X  2H 

SE 

(0.2366) 

1.195 

i 

1.92  X  2H 

BA 

(0.2946) 

1.942 

1.18  X  2H 

Isotropic 

(0.500) 

2.106 

1.09  X  2H 

Table  3.4:  Exact  characteristic  decay  lengths:  d  = 


Chapter  4 


THE  ANISOTROPIC  CASE 


The  techniques  described  in  the  previous  chapter  will  now  be  extended  to  the  general¬ 
ized  fourth-order  problem  for  the  anisotropic  strip.  First,  for  the  reader’s  convenience, 
we  will  recall  several  results  that  were  derived  in  Chapter  2.  Several  estimates  will 
then  be  presented  and  the  results  discussed. 

4.1  Formulation 

From  (2.21),  the  nondimensional  governing  equation  for  the  anisotropic  strip  is  given 

by 

2  12 

~  d"  ~  ~^,ivvn  d"  ~  ('^•^) 

^2  Cl  C3 

From  (2.34)  and  (2.35),  the  final  reduced  form  of  the  first  and  second  Conservation 
Laws  are 


62 
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J  =  I  [<!>%  +  d-q  (4.2) 

Liz  Cl 


and 


J/*  ^  4 

[2^, »>»?»»  +  +  —<l>,((v<f>,(vv]  dq,  (4.3) 


respectively.  From  (2.61),  the  nondimensional  physical  strain-energy  is  given  by 

2^12  /■! 

From  (2.65)  and  (2.72),  the  energy  functionals  Ei{z)  and  Eiiz)  are  given  by 

^  Ir}-^^'^^  ~  +  ^f,„]  (4.5) 

^  ~  ~  (4.6) 


From  (2.69)  and  (2.73),  the  line  integral  representations  for  these  energies  are 

Jf  12 

'  [-<t><f>,m  +  dq,  (4.7) 

i>x  Cl  62 

E2{z)  =  -  dq- 

JLx 


(4.8) 
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Furthermore,  it  is  assumed  that 


Cl 


1 


^  ,2  -2  ’ 
t2 


4 


(4.9) 


so  that  the  integrands  in  the  energies  Ei{z)  and  £^2(2)  are  positive-definite  quadratic 
forms  in  their  arguments. 


4.2  Basic  Energy  Estimate 

For  the  basic  energy  estimate  only  one  of  the  energies  Ei(z)  is  used.  As  in  (3.18)  for 
the  orthotropic  case,  let  the  function  F{z)  be  defined  by 

fOO 

F{z)  =  Eiiz)  +  2k  J  Ei{s)  ds,  (4.10) 

where  A:  >  0  is,  as  yet,  an  unspecified  constant.  From  (4.10),  we  have  that 

F'{z)  -1-  2kF{z)  =  E[{z)  +  4k^  Ei(s)ds.  (4.11) 

Returning  to  (3.20)  and  letting  /  = 

yields  an  expression  for  E^{z)  in  terms  of  a  line  integrzd  as 

=  -  J.  +  <!>%„]  (4.12) 

JLg  t2  Cl  C3 

Similarly,  letting  /  =  (3-20)  yields 

~  77  L  (4-13) 

dz  jRi  Cl  C2 
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which  when  integrated  using  the  divergence  theorem  and  boundary  conditions  gives 
an  expression  for  J  Ei{s)ds  in  terms  of  a  line  integral  as 

poo  t  12 

/  Ei{s)ds  =  - 

Jz  JRx  ^1  ^2 

=  -  W.£f  “ 


iV- 

Lti  t3 


(4.14) 


Returning  to  (4.11),  we  have  that 


F\z)  +  2kF{z)  =  -  /  -  ^<l>,ii<t>,in 


63  61 

.,2  42  2jfc2  -  8A:2  ,  ,  ^ 

-  4A:  (f>\ - —<f>%  +  dT) 

ti  t2 


-  Ji^)- 


(4.15) 


As  before,  we  want  to  find  values  of  k  such  that  J{z)  <  0.  The  steps  that  follow  are 

modelled  essentially  after  Horgan  [4],  except  that  here  we  use  the  energy  functional 

(4.5)  (rather  than  the  physical  energy  used  in  [4])  and  a  sharper  Wirtinger  inequality 

than  that  employed  in  [4]. 

If  we  define  Ei(z)  =  /  2WidA^  where 

J  Rx 

=  't’%  +  -  UiVt’Ar,  -  .  (^.16) 

£2  £3 

then  from  (4.15)  we  obtain 

.l(^)  =  -  pw,  -  41^(4  -  <‘1-  (4.17) 

We  can  write  2Wi  as  a  quadratic  form  as  s^Bs,  where 
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^,vv 


and  B  = 


(4.18) 


1  1 

Vita  2ei 


By  virtue  of  (4.9),  we  know  that  2Wi  is  positive-definite  and  we  have  the  result  that 

m  >  A„  +  2^;^,),  (4.19) 

where  A,„  >  0  is  the  minimum  eigenvalue  of  B  (see  [27]).  This  leads  to  the  inequality 

Jiz)  <  -  I 

J  Zjx 

-  4t=(^;'  -  dv  (4.20) 

Zti  t2 

=  -  L  +  2A„^f5, 

JLx  dAfn 

2/2  j2  .  ,  ,  4A:^  2l  j  /A  n^\ 

—  +  - 4>,v‘f>X  ~  T~^  ]  (4*21) 

Cl  €2 

<  Ml,  +  2A„,^i,  - 

8P  .  .  4k^  ,  _ 


SP  ,  ,  4k^  ,21  J 

+  -  T“^  ] 

€2 


(4.22) 


We  now  use  the  following  two  Wirtinger  inequalities  obtained  by  choosing  to  =  ^  in 
(3.12)  and  w  =  in  (3.11), 


—  id?  dn. 

(4.23) 

(4.24) 

respectively.  We  shall  set  h  =  2  (i.e.  equal  to  the  strip  width)  later.  Thus  we  obtain 
from  (4.22)  that 
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H - -  -T— ^  ]  drf.  (4.25) 

Observe  that  the  first  three  terms  in  the  integrand  may  be  written  as  a  quadratic 
form  as  Q  =  d'^Dd,  where 


UJ> 

1 _ 

and  D  = 

(2A^  -  4t^)  f 

= 

(4.26) 


Here,  D  is  a  symmetric  matrix  and  for  convenience  the  subscript  m  has  been  dropped 
from  the  eigenvalue  A.  If  we  choose  the  constant  k  such  that  the  eigenvalues  of  D  are 
positive,  i.e.  D  is  positive-definite,  then  we  obtain 

Q  >  a;(d'^d)  =  u}{<j>\  -1-  <l>\)  (4.27) 

where  a;  >  0  is  the  minimum  eigenvalue  of  D.  Equation  (4.25)  then  leads  to  the 
inequality 

r  Ah* 

j{z)  ^  K^fe  +  ^f.,) - dv-  (4.28) 

Dropping  the  first  term  and  using  the  inequality 

which  is  obtained  from  (3.11)  with  the  choice  w  =  we  then  have 


(4.30) 
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We  now  choose  k  to  be  the  positive  root  of  the  equation 


4k^ 


=  0 


(4.31) 


which,  since  A  >  0,  is  consistent  with  the  earlier  assumption  that  iw  >  0.  This  is 
the  largest  positive  value  for  k  such  the  right-hand  side  of  (4.30)  is  nonpositive. 

The  problem  may  now  be  simplified  using  the  following  notation.  Let 


k^  = 


Att* 


a  = 


2ei 


(4.32) 


The  matrix  D  in  (4.26)  can  then  be  written  as 


D  = 


Att* 


1-s  f 


-2-  2 -as 

<3 


(4.33) 


where  h  =  2  has  been  used.  The  new  notation  now  implies  that  in  order  to  solve 
for  the  positive  root  k  of  (4.31),  we  must  solve  for  the  positive  root  s  of  the  equation 


s^  =  2i/ 


(4.34) 


where 


is  the  minimum  eigenvalue  of 


2w 


(4.35) 


E  = 


1-s  f 
€2 


-2-  2 -as 

€2 


(4.36) 


Solving  for  i/  in  the  equation 


det\E  —  vl\  =  0, 


(4.37) 
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and  returning  to  (4.34),  we  seek  the  positive  root  of 


5"  =  3  -  (a+l)s  -  ^{(a-l)s  -  1}2  +  .  (4.38) 

This  is  equivalent  to  seeking  the  smallest  positive  root  of  the  quartic  equation 

s*  +  2(a  +  l)s^  +  (4a  -  6  -  4)^^  “  4(a  +  2)s  +  8  =  0,  (4.39) 

^2 

where,  from  (4.32),  a  =  Returning  to  (4.32)  thus  yields  the  following  value  for 
the  nondimensional  decay  rate  k, 


k  =  (4.40) 

where  s  is  the  smallest  positive  root  of  (4.39)  and  A  is  the  minimum  eigenvalue  of 
(4.18),  From  (4.18),  it  can  be  shown  that  A  is  the  smallest  root  of  the  characteristic 
polynomial 


2el 


+  =  (4-41) 


2e^ 


1_ 

2el 


With  (4.40)  as  our  value  for  k,  (4.15)  implies 


F{z)  <  F{0)e-^^\  z  >0.  (4.42) 

Following  exactly  the  arguments  given  in  Chapter  3  (in  (3.42)  -  (3.47)),  with  Ei{z) 
now  representing  the  anisotropic  energy  functional  (4.5),  we  obtain 

Ei{z)  <  2Ei{0)e-^’‘%  z  >  0.  (4.43) 

Again,  this  can  be  shown  to  lead  to  pointwise  estimates  for  the  stresses  r  such  that 


!’■(«.  >?)l  <  Ke-l* 


for  fixed  t). 


(4.44) 
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In  the  original  variables,  we  obtain  a  decay  estimate  for  r  with  exponential  decay 
rate  where 

F  =  (4.45) 

While  this  estimate  does  not  produce  an  explicit  formula  for  the  decay  estimate  in 
terms  of  the  parameters  ei,  €2,  £3,  the  calculation  of  A:  is  a  straightforward  procedure, 
involving  the  solution  of  the  quartic  and  cubic  algebraic  equations  (4.39)  and  (4.41) 
respectively.  Numerical  results  are  discussed  at  the  end  of  this  chapter. 

This  energy  method  may  also  be  carried  out  using  the  nondimensional  physical 
energy  Ep(z)  given  by  (4.4)  instead  of  the  mathematical  energy  Ei{z),  (Note  that  the 
analysis  in  [4]  maJces  use  of  the  physical  strain  energy.)  The  line  integral  representation 
of  Ep{z)  is  easily  obtained  by  using  (4.7)  and  observing  that 


(4.46) 


The  steps  that  follow  using  Ep{z)  are  essentially  the  same  as  those  presented  in  the 
preceding  analysis  for  E\{z),  the  main  dilference  being  that  the  problem  now  involves 
the  constants  ySp,  in  addition  to  the  parameters  cj,  €3  and  €3.  The  end  result  is  that 


\^22/  2y/2  V^22/ 


(4.47) 


where  s  satisfies  the  quartic  equation  given  by  (4.39)  but  now  a  =  — ,  and  A 

^yPiip22 

is  the  smallest  root  of  the  characteristic  polynomial 


^  ~y/^2 


+  -^u 


_  f  f  /^66  ^12 

V ^11^22  2  C2C3 


^12^66  X  _  _1_  _  _l_i  _  Q  4g\ 

2^11^22^  2el  2er  '  ^  ’ 
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Numerical  results  for  this  Ep{z)  formulation  were  calculated  and  compared  with 
the  results  jErom  the  Ei{z)  formulation  (see  Discussion  of  Results  section).  For  all 
of  the  materials  considered,  the  difference  between  the  two  results  is  rather  insignifi¬ 
cant.  The  mathematical  energy  Ei{z)  involving  the  parameters  ei,  62  and  €3  only  is 
algebraically  simpler  to  use  than  the  physical  energy  Ep(z).  An  advantage,  however, 
to  using  the  energy  Ep(z)  is  that  the  results  will  be  applicable  to  all  physical  mate¬ 
rials,  while  the  results  using  Ei{z)  depend  upon  the  condition  (4.9)  being  satisfied. 
This  condition,  which  guarantees  the  positive-definiteness  of  the  mathematical  energy 
norm  Ei{z),  fails  to  hold  for  some  materials  (see  Discussion  of  Results  section). 


4.3  Nonlinear  Optimization  Estimate 

This  estimate  makes  use  of  three  weighted  arithmetic-geometric  mean  inequalities, 
introducing  three  new  parameters  <x,  ^  and  7  into  the  problem.  Using  the  same 
function  F{z),  given  by  (4.10),  as  in  the  basic  energy  estimate,  we  recall  from  (4.15) 
that 

=  -J  [<f>%  + 

•'Lfx  C2 

t3  ti 

-  d-n-  (4.49) 

ti  t2 


We  seek  positive  values  for  k  such  that  J{z)  <  0.  Prom  the  arithmetic-geometric 
mean  inequalities  (3.15)  and  (3.16),  we  have  the  following  three  results 
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<  |^I^,J,II^.„I  <  +  J't’l,),  (4.51) 

-  <  ^(t€,  +  ^'fx)’  (4.52) 


where  a,  P  and  7  are  positive  weights  to  be  determined.  Using  (4.50)  -  (4.52)  in 
(4.49)  we  obtain 


Completing  squares  in  (4.53)  leads  to  the  inequality 

-  -L  l(‘“R){^«+(r^^  “  0%^' 

“ (^R + 0  (r  ^ 

which,  on  discarding  the  first  term,  yields 

-  ~L  I(‘“r)’‘’’"  ■*■  t  “  SR  “  ;9r) 

-  J^/\ 

\t2\) 


(4.55) 
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provided 

(l  -  >  0.  (4.56) 

Using  the  inequality  (4.24)  in  (4.55)  yields 


J(z)  < 


< 


provided 


1 

a|c2| 


-«b)  - 


drj  (4.57) 


>  0.  (4.59) 


Using  the  Wirtinger  inequalities  (4.23)  and  (4.29),  we  get  from  (4.58)  that 


(4.60) 


provided 


(4.61) 


We  now  seek  the  maximum  positive  k  such  that  J{z)  <  0  i.e.,  from  (4.60),  the 
maximum  positive  k  such  that 


(4.62) 
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It  is  clear  that  the  maximum  value  for  k  occurs  when  equality  is  chosen  in  (4.62), 
which  upon  using  the  quadratic  formula  yields 


+ 


(-r)( 


|e2|  2e^) 


+  41 


Equation  (4.63)  guarantees  a  positive  value  for  due  to  the  conditions  (4.56)  and 
(4.61),  which  in  turn  guarantees  a  real  value  for  k.  We  observe  that  the  formula 
for  k  given  by  (4.63)  involves  the  parameters  a,  ^  and  7  which  have  not  yet  been 
determined.  It  remains  to  choose  positive  values  for  these  parameters  such  that  k  is 
maximized  and  such  that  the  conditions  (4.56),  (4.59)  and  (4.61)  hold.  Since  k  given 
by  (4.63)  is  a  nonlinear  function  of  the  parameters  a,  ^  and  7,  and  the  constraints 
(4.56),  (4.59)  and  (4.61)  also  involve  nonlinear  functions  of  these  parameters,  this 
formulation  is  called  the  nonlinear  optimization  estimate. 

We  observe  that  the  constraints  (4.56)  and  (4.61)  may  always  be  satisfied  by  a 
proper  choice  of  a  and  The  constraint  (4.59),  however,  involves  k  as  well  as  a,  /?, 
7,  and  may  or  may  not  be  satisfied  by  the  choice  of  k^  in  (4.63).  If  a,  ^  and  7  are 
chosen  such  that  the  constraints  (4.56),  (4.59)  and  (4.61)  are  satisfied  with  k  given 
by  (4.63),  then  this  value  of  A:  is  a  valid  estimated  decay  rate.  If  (4.56)  and  (4.61)  are 
satisfied  by  this  choice  of  k  but  (4.59)  is  not,  then  instead  of  (4.63),  we  choose 


(ei  o|62|  ^Ical)  ((:5fe 


TT* 


+  1)]  4/l2’ 


(4.64) 
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and  so  equality  holds  in  (4.59).  It  is  easily  verified  that  k  given  by  (4.64)  satisfies 
(4.62), 

Equations  (4.63)  and  (4.64)  give  two  separate  formulae  for  the  estimated  decay 
rate  dependent  on  the  choices  of  a,  p  and  7  and  valid  under  differing  conditions. 
Furthermore,  while  there  are  many  choices  for  a,  ^  and  7  which  yield  a  valid  estimated 
decay  rate,  we  need  a  method  to  determine  the  proper  choices  of  a,  ^  and  7  which  yield 
the  optimal  decay  rate.  The  two  formulae  (4.63)  and  (4.64)  lead  to  two  optimization 
formulations  which  will  now  be  analyzed  to  determine  the  largest  value  for  k. 

4.3.1  Nonlinear  Optimization  Analysis 

Of  the  three  constraints  (4.56),  (4.59),  and  (4.61),  that  must  be  satisfied  for  the  non¬ 
linear  optimization  estimate,  two  of  them,  (4.56)  and  (4.61),  may  always  be  satisfied 
by  proper  choices  of  a  and  The  remaining  constraint  (4.59)  may  be  equivalently 
expressed  as  the  following  two  reduced  conditions. 


/9k3l)  -  ®’ 

-  —)  ^ . 


(4.65) 

(4.66) 


where  (4.65)  is  necessary  in  order  for  k^  to  be  nonnegative.  Prom  (4.56)  and  (4.61), 


a  <  Ical 

1 

-  > 
a 

1 

hi’ 

(4.67) 

P  <  ksl 

1 

ksl’ 

(4.68) 

we  have 
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which  in  turn  yield 


(ei  aleal  ^jeal)  ^  (ci  ej  e§)  ’ 


(4.69) 


where  we  observe  that  right-hand  side  of  (4.69)  is  positive  due  to  the  positive-definite 
assumption  (4.9)  on  the  the  energy  functional  Ei{z).  As  long  as  the  right-hand  side 
of  (4.69)  is  strictly  positive,  we  are  guaranteed  that  there  is  some  choice  of  a  and 
for  which  (4.65)  is  satisfied.  The  remaining  constraint  (4.66)  may  or  may  not  be 
satisfied  by  a  proper  choice  of  7,  and  thus  determines  which  of  the  two  formulae  (4.63) 
or  (4.64)  is  valid.  We  then  have  the  following  two  optimization  formulations: 


Optimization  Problem  1: 


Maximize  F(a,  /3,  7)  =  “  j^)  (]^  +  ^) 


subject  to 


and 


TT^ 

(-h)'i 

(hr 

-V 

2€,) 

+  4^1 

Nl) 

r) 

(■■ 

a  ) 

■  1^21) 

1  >  » 

(4.71) 

■  ksU 

1  -  ° 

(4.72) 

( - 

\ei  a 

1 

k2|  " 

1  ) 

1  ^  0 

(4.73) 

/'JL+ n 

\(-- 

1 

1  ^ 

(  7hl  \ 

Vz\) 

VI'jI 

«l«2| 

VI  +7hl/ 

1 

1  V/" 

2 

—  A  1 

(  ^ 

1 _ _ 

"l  (l  - 

>  (A  JA^ 

a|c2| 

\1+7|C2|/ 

1^21 

jv  n; 
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Optimization  Problem  2: 

Maximize  G(u,  ft  t,)  =  ^  ^ 

subject  to 


7|c2| 


+  71^21 


) 


and 


(-h)( 


I  €2 1  2ei 


+  • 


1  _ 1  V/'  7N 

\ei  a 


(r 


(4.75) 


-a 

>  0 

(4.76) 

-a) 

>  0 

(4.77) 

>  0 

(4.78) 

1 

1  W  7N  \ 

"|C2| 

^Icsl/  VI  +7^21/ 

-V  -  4  fl 

<  0.  (4.79) 

>\)  1 

Condition  (4.74)  results  from  using  (4.63)  in  (4.66)  and  condition  (4.79)  is  the  in¬ 
equality  that  results  when  (4.74)  is  not  satisfied. 

To  maximize  the  functions  F{a,  )9,  7)  and  G{a,  7),  we  now  look  at  their 
partial  derivatives  with  respect  to  the  parameters  a,  ^  and  7.  Prom  (4.70),  it  can  be 
shown  that 

F^a  <  0,  <  0,  F,y  <  0.  (4.80) 

Similarly  from  (4.75),  it  can  be  shown  that 


Ga  >  0. 


G^a  >  0, 


G,p  >  0, 


(4.81) 
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Thus  to  maximize  E(o!,  fl,  7)  we  seek  the  smallest  possible  values  of  a,  P  and  7  while 
to  maximize  G{a,  yS,  7)  we  seek  the  largest  possible  values  of  a,  ^  and  7.  It  will  now 
be  shown  that  in  fact  the  choices  of  a,  /5,  7,  call  them  o:*,  and  7*,  which  maximize 
F  also  maximize  G. 

For  fixed  values  of  a  and  ^  such  that  (4.71)  -  (4.73)  are  satisfied,  F  is  a  decreasing 
function  of  7  and  hence  the  choice  yp  which  maximizes  F  is  the  smallest  positive 
7  such  that  (4.74)  is  still  satisfied.  For  fixed  a  and  P,  the  smallest  positive  7  for 
which  (4.74)  is  satisfied  occurs  when  equality  is  taken,  yielding  a  cubic  polynomial  in 
7  with  one  positive  root  ')p.  At  this  value  of  7  =  7'p,  the  functions  F(a,  P,  7)  and 
G{a,  /?,  7)  are  identical  (recall  that  (4.74)  results  from  (4.66))  and  both  are  valid 
decay  estimates.  As  7  decreases  below  7]^,  (4.79)  is  satisfied  and  G  becomes  the  valid 
decay  estimate.  With  a  and  yS  still  fixed  such  that  (4.76)  -  (4.78)  are  satisfied,  G  is  an 
increasing  function  of  7  and  the  choice  7^  which  maximizes  G  is  the  largest  positive 
7  such  that  (4.79)  is  satisfied,  i.e. 

=  If  =  7*.  (4.82) 

Furthermore,  we  have 

f  (a,  7")  s  (?(«,  7.).  (4.83) 

So  for  fixed  values  of  a  and  the  optimal  choice  7*,  which  maximizes  both  F  and 
G,  is  foimd  from  computing  the  positive  root  of  a  cubic  polynomial  and,  since  F  and 
G  are  identical  at  this  choice  of  7*,  k  is  obtained  from  either  (4.70)  or  (4.75),  using 
h  =  2.  It  then  follows  that  over  the  full  range  of  a,  and  7,  the  choices  a*,  and 
7*  which  maximize  F  also  maximize  G  and  that  the  two  functions  are  identical  at 
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these  values. 

While  (4.70)  aud  (4.75)  yield  two  explicit  formulae  for  a  decay  estimate,  they  are 
in  terms  of  three  parameters  that  must  be  optimally  chosen  and  a  numerical  procedure 
was  necessary  to  do  this.  A  program  was  implemented  to  run  through  the  possible 
ranges  of  the  parameters  a  and  /S,  obtained  from  (4.71)  -  (4.73)  or  equivalently  (4.76) 
-  (4.78),  and  compute  an  optimal  value  k  for  each  pair  of  a  and  (in  the  manner 
previously  described),  and  thus  an  optimal  k  over  the  entire  range  of  a,  ^  and  7.  The 
optimal  decay  estimate  k  in  the  nondimensional  domain  was  then  transformed  back 
into  the  original  domain  so  that 

=  (4,84) 

and  the  stresses  t  have  an  exponential  decay  rate  of  These  results  axe  shown  at 
the  end  of  the  chapter.  It  is  observed  that  while  this  nonlinear  optimization  estimate 
produces  explicit  formulae  for  the  estimated  decay  rate  in  terms  of  parameters  which 
may  be  optimally  chosen  numerically,  the  formulae  (4.70)  and  (4.75)  will  be  different 
for  each  material  since  the  parameters  o,  an  7  must  be  optimally  chosen  in  each  case. 
This  is  in  contrast  to  the  formulae  (3.51),  (3.55)  and  (3.102)  obtained  in  Chapter  3  for 
orthotropic  materials.  These  formulae  did  not  change  from  one  material  to  another. 
Ideally,  we  would  like  to  determine  the  optimal  choices  of  a,  and  7  analytically  in 
terms  of  the  material  parameters  e\,  €2  and  C3.  Then  the  formulae  (4.70)  and  (4.75) 
would  be  explicit  in  terms  of  Ci,  €2  and  €3  and  valid  for  all  materials  satisfying  (4.9). 

If  we  now  consider  the  orthotropic  limit  for  which  ^  “*■  0,  this  nonlinear 
optimization  formulation  reduces  to  the  following: 
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Problem  1 


= 


’’•Ml 

+  :77ta/t^  +  4 


2h'^  V  iej 


where 


—  >  0  and  —  4  >  0. 

ei  -  4ei 


(4.85) 


(4.86) 


k  = 


TT 

2hy^^ 


(^1  +  let?  - 


(4.87) 


Problem  8 


where 


1  3 

—  >0  and  —X  —  4  <  0. 
ei  4ef 


(4.88) 


(4.89) 


=>►  k 


n  ^/3 

2hy/€l'  ^  4  ■ 


(4.90) 


On  comparing  (4.87)  amd  (4.90)  with  (3.49),  we  see  that  in  the  orthotropic  limit  the 
nonlinear  optimization  estimate  reduces  to  the  analytic  estimate  (3.49)  developed  in 
Chapter  3.  This  is  not  unexpected  as  the  methods  used  to  develop  the  two  estimates 


are  similar. 
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An  alternative  formulation  for  this  nonlinear  optimization  method  may  be  ob¬ 
tained  using  the  physical  energy  Ep(z),  given  by  (4.4),  instead  of  the  energy  func¬ 
tional  Ei(z).  This  leads  to  a  nonlinear  optimization  estimate  in  /our  parameters  a, 
/3,  7  and  S  characterized  by  the  following  two  problems: 


Optimization  Problem  1: 


\^12\S 


F(a,  (l  -  ^ 

Yi  - — - 

|C2|  V  ^11^22  J  \IC2I  2  •y ^11  ^22/ 


^66  \ 


TT 

2h^ 

+  4 


subject  to 


V  1^2!  V^11^22/  V 


12 


ksl  ^11^22 


)]' 


(4.91) 


OtThd 


/,  _  C.  \MS  \ 

V  k2|  V ^11022/ 

(l  _  Jt _ 1^2!  \ 

\  l^sl  ^•n/^1i/?22/ 

(gee  _  J _ l_\ 

yj ^11^22  ^\^2\ 


1^12!^  \  ( _7_  ,  ^66  \  /  ^66 

V^ll/322/  \k2|  2-v/^11^22/  \VAi^22 


>  0 
>  0 
>  0 


1 

aleal 


(4.92) 

(4.93) 

(4.94) 


M  /  7|e2| 

^N/  \H-7hl 


) 


,  1  /  ^ee _ ^ 

4  \ V^11^22  “h2| 


-  4 


y/^2)  I 


1  Vf  7n  y 

Vl  +  7hl/ 

1 _ ^■^=]  >  0, 

jcal  by/j^n^J 


(4.95) 
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Optimization  Problem  2: 


Maximize  G(a,  7,  5)  = 


/  ^66 

4h^  \\/Pn^22 


1 

a|e2| 


M  /  7|e2|  \ 
M)  \1  + 71^21/ 

(4.96) 


subject  to 


(■ 

a  \ 

1^2!  y/ 011022) 

> 

0 

(4.97) 

?  IAjI  'i 

> 

0 

(4.98) 

V 

jcsl  by/ ^11022/ 

(  006 

1  1 

j  > 

0 

(4.99) 

\V  011022 

oIcjI  0\t3\) 

(  J-j.  (  ^66  _  _1 _ [  7^2! 

\  |€2|  V^n^)  \\^2\  •2-v/;0ii^22/  VV/^11^22  “I^^I  /^ksl/  +  7^21/ 


1  /  1^66 _ 1 _ 1  v/  7^2!  y 

^4  \V ^11^22  oc\€2\  ^leal/  \H-7|e2l/ 

\M  \  <  „ 

V  1^2!  V 011^22)  \  Nsl  011022/ 


(4.100) 


Numerical  results  were  calculated  for  this  Ep{z)  formulation  and  compared  with  the 
Ei{z)  formulation.  As  was  the  case  previously  for  the  basic  energy  estimate,  the 
difference  between  the  two  energy  formulations  was  insignificant  and  it  is  clear  that 
the  Ei(z)  formulation  is  much  simpler  to  work  with,  keeping  in  mind  however  that 
the  physical  energy  formulation  has  a  wider  range  of  applicability. 
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4.4  Higher- Order  Energy  Estimates 

In  this  section  we  adapt  the  conservation  law  approach,  developed  for  orthotropic 
materials  in  Chapter  3,  to  anisotropic  materials,  making  use  of  the  higher-order  energy 
functional  E2{z).  Thus  for  the  estimates  that  follow,  both  of  the  anisotropic  energies 
Ei{z)  and  E2{z)  given  by  (4.5)  and  (4.6)  are  used.  Let  the  function  F(z)  be  defined 
by 

F{z)  =  [E2{z)  +  mEi{z)]  +  j  [-^2(5)  +  mEi(s)]  ds,  (4.101) 

where  m  and  k  are  both  positive  constants  to  be  determined.  With  F{z)  defined  in 
(4.101),  we  have 

F'{z)  +  2kF{z)  =  E'^{z)  +  mE[{z)  +  4k^  [i;2(5)  +  m£i(s)]  ds.  (4.102) 

t 

We  now  wish  to  express  the  quantities  E2{z)  and  J  E2(s)ds  as  line  integrals,  as  was 

roo 

done  previously  for  E'i{z)  and  J  ^i(s)ds  in  the  basic  energy  estimate  of  section  4.2. 

From  its  definition  in  (4.6),  and  using  (3.20),  we  obtain  a  line  integral  represen¬ 
tation  for  E2{z)  as 


+  —^%r,  - 

€2  €1 


‘‘I-  (‘•“'3) 

C3 


Recalling  (4.8),  we  obtain 


^2(2)  =  ^  ^ 


(4.104) 
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which,  when  integrated  and  combined  with  the  divergence  theorem  and  boundary 
conditions,  yields 


(4.105) 


We  observe  that  (4.105)  is  identical  to  (3.60)  obtained  for  the  orthotropic  case  in 
section  3.4,  i.e.  the  quantity  /  E2{s)ds  has  the  same  line  integral  representation 
regardless  of  whether  ■£'2(5)  is  the  ortftotropic  functional  (3.8)  or  the  anisotropic  func¬ 
tional  (4.6).  This  is  because  the  line  integral  representation  (4.8)  for  the  anisotropic 
energy  functional  E2{z)  is  independent  of  any  material  parameters.  Substituting 
(4.12),  (4.14),  (4.103)  and  (4.105)  into  (4.102)  results  in 


F\z)  -I-  2kF{z)  =  -  f  -b 

JLx  ^2  ^3 


2\42 


+  +  (m  -  +  (m  -  2k^)<l> 

2»71  ,  ,  2m  ,  ,  m  ,2  .,2/2 

€2  C3  €1 

.t2  I  /  8Pm  ,  ,  2Pm  ,2,  , 

+  4A:  -b  — - — <f>,((f>,r, - ; — (pj  dr] 

€2 


(4.106) 


We  will  now  investigate  several  methods  to  obtain  positive  values  for  m  and  k  such 
that  J{z)  is  nonnegative.  This  in  turn  yields  the  exponential  decay  estimate 

F{z)  <  F(0)e-2*%  z  >  0,  (4.107) 


which,  by  following  exactly  the  arguments  outlined  in  (3.85)  -  (3.89),  leads  to  bounds 
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on  the  energy  functionals  of  the  form 


Ei(z)  <  2  [Ei(0)  +  -^2(0) 

L  m  J 


-2kz 


z  >  0, 


(4.108) 


Eiiz)  <  2  [£2(0)  +  mEi{0)]e-^'‘%  z  >  0, 


(4.109) 


where  ^'i(^)  and  E2{z)  now  represent  the  anisotropic  energy  functionals.  Again,  the 
results  (4.108)  and  (4.109)  can  be  used  to  derive  pointwise  estimates  for  the  stresses 
T  such  that 

|r(^,  7/)|  <  Ke~^^  for  fixed  rj.  (4.110) 

4.4.1  Quadratic  Estimate 


This  estimate,  while  modelled  after  the  conservation  law  approach  of  Chapter  3, 
does  not  actually  involve  either  of  the  anisotropic  conservation  properties  (4.2)  or 
(4.3).  It  does  however  use  the  combination  of  energy  functionals  Ei(z)  and  E2{z) 
that  was  used  in  Chapter  3  for  the  orthotropic  conservation  law  estimate.  The  result 
is  a  somewhat  simpler  quadratic  polynomial  for  m  in  contrast  to  the  the  sixth-order 
polynomial  developed  in  the  orthotropic  conservation  law  estimate. 

Starting  with  (4.106)  and  completing  squares,  we  obtain 


j(‘)  =  L 

JLg  €2  ^2  63  €3 

+  +  (>^  - 

,  ,  2m  ,  ,  .  ,2  ai2  i2 

-  +  —r,in  “  rn<f>,i 

t2  C3  Cl 

.,2  ,  ,  8Pm  ,  ,  2Pm  .3 ,  . 

+  4A:  - - — 4>W  dr) 


£2 


(4.111) 
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2 


-  X.  [  ~  i  ■  5)  ~ 

2  2  1 

€2  C3 

-  dt}.  (4.112) 

Completing  squares  again  yields 

U2_^2  ,  -  — *f,l  rf,  (4.113) 


—  4k^m(f>^^  +  4k^m(j>(f>^^(^  + 


C2 


2\J'2 


)4> 


which  upon  dropping  the  completed  square  term  gives 

+  m  (-  -  -  ^<l>.nn<f>,^v  “ 

\Cl  ^2/  ^3 

+  +  ^^^<f>,(<t>,r,  -  -7— drj.  (4.114) 

£2  ^1  J 

Next,  we  use  several  Wirtinger  inequalities  to  simplify  the  expression  for  J{z)‘  Using 
(4.23)  and  (4.24)  in  addition  to  (3.70)  from  Chapter  3,  we  obtain 


+{-(j-3)-“l« 


2 

m 


,iv 

'•  J 

.,1  II  SPm  .  ,  2m  ,  , 

+  4fc'‘m#,(^  +  - - <t>,m9iiv 

€-2  63 


dt]. 


(4.115) 


Using  the  arithmetic-geometric  mean  inequalities  (3.15)  and  (3.16),  we  have,  for  ar¬ 
bitrary  constants  7,  o:  >  0, 


CHAPTER  4.  THE  ANISOTROPIC  CASE 


87 


I  4k^m<l)<f>^^^  dr]  >  —  I  2k^m[a<j)^  +  —<t>%]  drj 
JL/z  JIjz  ^ 


where  the  last  inequality  follows  on  using  (3.14),  with  w  =  <f>.  Substituting  these 
inequalities  into  (4.115)  yields 


,2  2Pm\  2 


>  L  [{£-4-3 

f  «.2  2k^mh‘^  2k^mh*a  m'y\ 

h 


{”■£-3 


jVV 


m 


4fc^m 

'Td" 


+  ]<;-)•  (4.119) 


Again  using  (4.23)  and  (4.24)  we  obtain 
+  |m- 


2  21:2m  \  2 


a 


} 


■t: 


2  2k^mh^  2k^mh‘*a  m'f  4k^mh^  \  ^ 

”  R  ”  r 


47r*€i 


f  / 1  1  ^  4A:2m^2  02  4k'^mh?P\  ,2  1  j  /,i  10m 


It  now  follows  that  J{z)  >  0  if  the  positive  constants  m,  k,  a,  ^  and  7  can 
be  chosen  so  that  the  coefficients  of  and  in  (4.120)  are  all  nonnegative. 
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Thus  we  seek  m,  k,  a,  ^  and  7  such  that 


A:*  < 


(L. 

1 

1] 

TT^ 

Vci 

~ 

'  4j 

1  ^ 

* _ ft _ / 

K-S 

1  1  1 
€1  4  kalT 

V  hi; 


(4.121) 


(4.122) 


mh?  _  amh*  ^  mh^ 


(4.123) 


In  order  for  k^  to  be  positive,  the  right-hand  sides  of  (4.121)  -  (4.123)  must  also 
be  positive.  Recalling  our  fundamental  assumption  (4.9)  on  the  Cj,  we  see  that  the 
right-hand  side  of  (4.121)  is  indeed  positive.  The  remaining  two  conditions  (4.122) 
and  (4.123)  imply  the  following  bounds  on  7, 

RTTTTY  <  ^  <  l-l- 

It  is  readily  verified  that  the  right-hand  sides  of  (4.121)  -  (4.123)  are  strictly  monotone 
functions  of  the  Vciriables  m,  a,  7  and  /3.  Thus  the  maximum  choice  for  k  satisfying 
(4.121)  -  (4.123)  results  from  equating  these  three  expressions.  Equating  (4.122)  and 
(4.123)  yields  the  following  expression  for  m, 


1_1 
Cl  c? 


/ 1  1  1  \ 
yci  el  |e3|7j 
1  ah^  I  1  ^  ^  f  ^ 

w;  V 
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while  equating  (4.121)  and  (4.122)  results  in  the  following  expression  for  a, 


It  can  be  verified  that 

13  >  0  ^  F(7,  13)  >  0  (4.127) 

for  7  in  the  range  given  by  (4.124),  so  that  positive  values  for  ^  and  m  yield  a  positive 
value  for  a.  Substituting  (4.126)  into  (4.125)  gives  the  quadratic  equation  in  m 

Am*  +  Bm  +  C7  =  0  (4.128) 


C 


/i*  [ei  el 


The  roots  of  (4.128)  may  be  found  explicitly  from  the  the  quadratic  formula  as 
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The  inequality  (4.135)  results  from  (4.131)  which  gives  a  quadratic  inequality  for  ^ 
while  the  inequality  (4.136)  results  from  (4.132).  Of  the  three  constraints  (4.134)  - 
(4.136)  that  must  hold  for  this  quadratic  estimate,  (4.134)  can  always  be  satisfied  by 
a  proper  choice  of  7  while  (4.135)  then  gives  a  lower  bound  on  the  choice  of  ^  in  terms 
of  7.  Furthermore,  it  can  be  observed  that  the  left-hand  side  of  (4.136)  grows  without 
bound  for  /3  large  and  increasing,  while  the  right-hand  side  decreases  with  respect  to 
/?,  so  that  (4.136)  is  always  satisfied  for  ^  sufficiently  large.  Thus,  we  are  guaranteed 
some  choice  for  7  and  ^  that  yields  a  valid  estimated  decay  rate  k.  While  (4.133) 
gives  an  explicit  formula  for  k,  numerical  procedures  are  necessary  to  determine  the 
parameters  ^  and  7  that  maximize  k. 

Without  implementing  any  numerical  techniques  to  maximize  k,  we  observe  from 


(4.122)  and  (4.124)  that 

where  the  right-hand  side  of  (4.137)  is  positive  by  virtue  of  (4.9).  Transforming  this 
result  back  to  the  original  domain  and  setting  h  =  2  yields  the  following  bound  on 
the  estimated  decay  rate  k* 

We  recall  that  the  stresses  r  have  a  decay  factor  of  The  inequality  (4.138)  yields 
an  upper  hound  for  the  estimated  decay  rate  solely  in  terms  of  the  elastic  parameters 
of  the  material,  i.e.  independent  of  7  and  and  will  be  used  later  to  assess  the  range 
of  validity  of  the  method  given  here.  It  will  be  shown  that  the  quadratic  estimate 
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is  limited  in  certain  ranges  of  the  material  parameters.  For  this  reason,  numerical 
optimization  techniques  were  not  pursued. 

We  now  return  to  (4.121)  -  (4.123)  and  consider  the  orthotropic  limit  as 
^  ^  0.  Equating  the  three  orthotropic  versions  of  (4.121)  -  (4.123)  leads  to  a 

quadratic  equation  in  m  that  is  simpler  than  in  the  anisotropic  case  and  leads  to  the 
following  orthotropic  quadratic  estimate: 

=  4^  - 

provided 

i  -  2€,  <  (4.140) 

which  implies 

d  >  iV  «  0.639  .  (4.141) 

The  inequality  (4.140)  is  a  condition  that  is  necessary  for  m  to  be  a  positive  real 
quantity.  Equation  (4.139)  is  identical  to  the  analytic  estimate  (3.49)2,  however  the 
range  of  €i  given  by  (4.141)  is  more  restrictive  than  the  range  ci  >  0.433  required 
for  validity  of  the  analytic  estimate. 

4.4.2  Quartic  Estimate 

This  estimate  starts  with  J{z)  given  by  (4.106)  and  makes  use  of  the  anisotropic 
conservation  property  (4.2)  to  simplify  terms.  This  leads  to  a  quartic  equation  for 
the  undetermined  constant  m,  in  contra.st  to  the  quadratic  equation  for  m  which 
resulted  in  the  the  previous  section. 

The  conservation  property  (4.2)  may  be  rewritten  using  the  arithmetic-geometric 
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mean  inequality  (3.15)  as 

L  -  '!>%  +  -  f  4,1  ‘‘’I  =  -f,  (^.«5  <'-) 

JL/z  Lc2  J  L/x 

<  I  +  df)  (4.142) 

J  Lx  L  ^ 

for  tj  >  0.  Thus  we  have  from  (4.142)  that 

~^Jl,  -  ~  Jl,  ~ 

(4.143) 

Returning  to  (4.106)  and  using  (4.143)  leads  to  the  inequality 

J{z)  >  -  ^J>Ar,r,4>,iin 

-  —Mi.  +  [T^-^yi.  -  {*'=  "•+— )4 

+  ik'‘mil4i^  +  4I  <*')•  (4.144) 

^2  J 


Completing  squares  in  the  higher-order  cross  terms  results  in 


from  which,  upon  dropping  positive  terms,  we  get 
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-  L  [(e.  (l-£)£i 

+  (m-2k^-j^<f>%  +  (m  -  2A:2  +  y) 

-  —Mav  +  -  (4fc  m  +  —  j 


.,2  ,,  SPm  .  .  2Pm,2l  , 

+  dri,  (4.146) 

€2  Cl 


provided 


(■-£) 


(4.147) 


On  using  (3.66)  from  Chapter  3,  we  have 


m  2Pmi//i^l  .2  2m  m 

m-2k  +y - (3£)4  ”  —%r,n<l>Av  + 


SPm 


4  -  - 


£3  ZCi 

2Pm  ,,1  , 


/^.o  wm\  ,2  SPm  ,  .  2Pm,2l  ,  /. -...ox 

-  (^4l:^m  +  —  j  +  —^(f>,^4>,r, - ('^•1^®) 

where  1/  >  0  is  an  arbitrary  parameter.  Now,  using  the  arithmetic-geometric  mean 
inequalities  (3.15)  and  (3.16),  we  obtain 

for  arbitrary  S,  pi  >  0  .  Using  these  inequalities  in  (4.148)  yields 


4A:*m  -t-  -—  -|- 


wm  4k^mS 


h-{ 


2k^m  4Pm\  ,2!  ,  irix 
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Upon  using  (4.23)  and  (4.24)  in  the  last  two  terms  of  (4.151),  we  obtain 


J(!)  > 


L  l{j 


(1  -  ^)4 


+ 


{m  m 

2ei  n\e3\ 


2k^muh^  iim  k'^mh?'  k^mh"^  1  ,, 

<P 


^3 


j 


,TiV 


(f)<  leal  2£,7r! 

ik^mh^  muh^  4f^mh^S\  1  .  ,, 

- 


TT'* 


27r2 


On  using  (3.70)  we  obtain 


m  > 


L  l(S(; 


a 


+ 


+ 


(1  - 

{-(I-h) -“■(-* 


muh^  mh^ 

+  ■: - :r  + 


"‘V)) 


4>: 


)“*  4€i7r2  26\t2\ 


(4.153) 


provided 


( *.  (1  -  £)4  4)  -  "• 


(4.154) 


It  then  follows  that  the  integral  J{z)  given  by  (4.153)  will  be  positive  if  we  can  choose 
positive  constants  w,  i/,  8,  fj,,  m  and  k  such  that  the  coefficients  of  and 

are  nonnegative  and  such  that  the  conditions  (4.147)  and  (4.154)  are  both  satisfied. 
Since  (4.9)  implies  that  £§  >  €1,  these  last  two  conditions  may  be  combined  together 
to  yield  the  single  constraint 


r^\  ^  eie| 

2u)  -  el{el-e^y 


(4.155) 


which  shall  be  used  henceforth. 
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Thus  the  problem  is  to  find  positive  constants  iv,  u,  6,  fi,  m  and  k  such  that 


< 


< 


< 


/_! _ 

\2ei  /fleal  2xV 

m  /3  II  \ 

?U~Mj 

mvh^  mh^  mh^  ’ 
4ei7r2  25|c2|7r2 

47r^  f  1  1  On 

•Fk-gFii-ar  » 


(4.156) 


(4.157) 


(4.158) 


2u:  -  ^  ~  elid  -  Cl)  = 


(4.159) 


For  convenience,  we  now  choose  equality  in  (4.159)  to  obtain  a  valid  decay  estimate 
although  this  may  or  may  not  yield  the  optimal  decay  estimate.  Upon  choosing 
equality  in  (4.159)  and  introducing  the  notation 


\A  A  1.2  2 

=  F">  "  =  F^’  *  =  F^  ’ 

where  h  =  2,  the  inequalities  (4.156)  -  (4.158)  then  read 


< 


J _ 1 _ M\ 

2ei  p\e3\  4/j 


< 


M  /3  _ 

2  U  lealj 


(4.160) 


(4.161) 


(4.162) 


< 


('+a) 


(4.163) 
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The  problem  is  now  to  find  positive  Af,  A,  and  K  such  that  (4.161)  -  (4.163)  hold, 
where  6  and  n  are  arbitrary  positive  parameters  yet  to  be  determined. 

First,  we  observe  that  in  order  for  to  be  positive,  the  right-hand  sides  of  (4.161) 
and  (4.162)  imply  the  following  restrictions. 


I _ 1 _ M 

2ei  /ileal  4/  ^ 

3|  .  ' 

A*  <  2^^3\- 


(4.164) 

(4.165) 


Since  the  right-hand  sides  of  (4.161)  -  (4.163)  are  monotone  functions  of  their  argu¬ 
ments,  the  optimal  value  for  K  will  be  obtained  by  equating  these  three  expressions. 
This  leads  to  a  fourth-order  equation  for  M  as 


1 


(I)  + r)  "  (r 


n  M 

1  +  J - 1" 

4ei 


M  \  /  J _ 1_ 

26\e2\)  \2ei  /ileal 


-a) 


0. 


(4.166) 


Once  /I,  subject  to  (4.165),  and  6  have  been  chosen,  equation  (4.166)  may  then  be 
solved  numerically  for  M.  Of  the  four  possible  roots,  M  must  be  chosen  such  that 
the  constraint  (4.164)  is  satisfied.  K  may  then  be  found  on  equality  in  (4.161)  and 
A  in  turn  may  be  found  from  (4.163)  with  equality.  Due  to  the  complexity  of  this 
formulation  and  the  number  of  parameters  involved,  numerical  calculations  were  not 


pursued. 
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One  can,  however,  obtain  a  bound  on  this  quartic  estimate  as  was  done  for  the 


previous  quadratic  estimate.  Returning  to  (4.161),  we  observe  that 

_1 _ 1_ 

<  ^  -  A 

4  8ei  6e5 


(4.167) 


where  the  condition  (4.165)  has  been  used  in  the  last  step.  Referring  back  to  the 
notation  in  (4.160)  and  transforming  (4.167)  back  into  the  original  xi  —  X2  domain 
yields  the  following  upper  bound  on  the  estimated  decay  rate  k*, 


r  <  lj±  -±(h.y 

2y  8ei  bcl  \^22/ 


(4.168) 


where  the  stresses  r  have  a  decay  factor  of  This  bound,  solely  in  terms  of  the 
elastic  parameters  of  the  material,  will  be  used  later  to  assess  the  range  of  validity 
of  the  present  analysis.  It  will  be  shown  that  the  quartic  estimate,  like  the  quadratic 
estimate,  is  limited  in  certain  ranges  of  the  materi2il  parameters. 

Again,  it  is  easier  to  obtain  results  from  this  method  in  the  orthotropic  limit,  i.e. 
0.  The  inequalities  (4.156)  -  (4.158)  then  reduce  to 


as  — ,  — 

ta’  £3 


< 


u> 
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Taking  the  limit  in  (4.159)  also  gives 
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From  (4.169)  we  see  that  A;  is  a  decreasing  function  of  ui  and  therefore  choosing 
equality  in  (4.172)  will  yield  the  optimal  decay  estimate.  Upon  using  the  notation  in 
(4.160),  the  remaining  conditions  (4.169)  -  (4.171)  become 


< 

< 


3M 


< 


(4.173) 

(4.174) 

(4.175) 


In  order  for  to  be  positive,  (4.173)  implies  that 


M  < 


d’ 


(4.176) 


Due  to  the  monotone  nature  of  the  right-hand  sides  of  (4.173)  -  (4.175),  the  maximum 
value  for  K  will  be  achieved  by  equating  these  three  expressions.  This  leads  to  the 
following  fourth-order  equation  for  M, 


(D*  a + + ii)  (f  - 


The  roots  of  (4.177)  may  then  be  used  to  obtain  K  by  tciking  equality  in  (4.173), 
where  it  is  observed  that  the  smallest  value  for  M  yields  the  largest  value  for  K. 
The  fourth-order  equation  (4.177)  was  solved  numerically  using  MACSYMA  for  the 
set  of  materials  shown  in  Table  3.1  from  Chapter  3.  For  each  of  the  materials,  the 
smallest  positive  root  M  gave  positive  values  for  A  and  and  hence  gave  the  desired 
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maximum  iif  satisfying  (4.173)  -  (4.175).  Using  (4.160)  and  transforming  this  result 
back  to  the  original  domain  yields 

where  the  stresses  r  have  a  decay  factor  of  Although  this  estimate  does  not 
yield  an  explicit  formula  for  the  decay  rate  in  terms  of  the  elastic  constants,  it  does 
yield  a  numerical  estimate.  These  results  are  shown  in  Table  4.1  together  with  the 
conservation  law  estimate  of  section  3.4.  It  is  seen  that  the  latter  method  yields  a 
better  estimate. 

4.5  Exact  Decay  Rates 

Exact  solutions  of  the  anisotropic  equation  (4.1)  subject  to  the  boundary  conditions 
(2.23)  -  (2.25)  are  much  more  complicated  and  difficult  to  obtain  than  the  solutions 
of  the  orthotropic  equation  (3.3)  subject  to  the  same  boundary  conditions.  As  in 
section  3.5,  we  seek  solutions  of  the  form 

(f>  =  7  =  constant^  (4.179) 

which  leads  to  the  following  eigenvalue  problem 

F""  +  -7F"'  +  +  -^^F'  +  7‘‘F  =  0,  on  (-1,1) 

€■2  €3 

F(-l)  =  F(l)  =  0,  (4.180) 

F'(-l)  =  i^'(l)  =  0, 

which  is  the  anisotropic  generalization  of  (3.104)  Seeking  solutions  to  (4.180)  of  the 
form 


F  = 


ft  constant, 


(4.181) 
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leads  to  the  characteristic  polynomial 

P{fi)  =  +  -fi  +  1  =  0  (4.182) 

Cj  Cl  C3 

whose  roots  fi  axe  complex  and  occur  in  conjugate  pairs  (Equation  (4.182)  is  actually 
the  nondimensional  version  of  (2.47)).  Let  the  roots  fi  of  (4.182)  be  denoted  by 

fil,2  =  Pi  ±  *91,  P3,4  =  P2  ±  *92,  (4.183) 

where  pa,  qa  depend  on  e,-  (i  =  1,2,3).  Substituting  the  general  form  for  F  given  in 
(4.181)  into  the  boundary  conditions  given  in  (4.180)  results  in  the  following  eigen- 
condition  for  7  (see  Choi  and  Horgan  [6]), 

[(Pi  -  P2)^  +  (91  +  92)^]  cos  2  {qi  -  92)7  - 

[(Pi  -  P2)^  +  (91  -  92)^]  cos  2{qi  +  92)7  =  49192  cosh  2  (pi  -  p2)'f.  (4.184) 

The  roots  of  (4.184)  are  complex  in  general  and  depend  on  the  three  eltistic  parameters 
€1,  £2  and  £3.  Recently,  in  [19],  this  eigencondition  has  been  simplified  using  a  spherical 
representation  to  yield  a  reduced  eigencondition  for  a  modified  eigenvalue  C  in  terms 
of  two  parameters  xp  and  u>  which  may  be  determined  from  the  elastic  constants  of  the 
material  (in  [19],  the  quantity  w  is  denoted  by  6).  The  actual  eigenvalue  7  is  related 
to  the  modified  eigenvalue  (  by  the  scaling 

C  =  P7,  (4.185) 

where  p  is  a  third  material  parameter  which  may  be  determined  from  the  elastic 
constants  of  the  material.  The  reduced  eigencondition  can  be  solved  for  various  values 
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of  xj)  and  u,  see  [19].  By  interpolating  from  tables  in  [19]  and  scaling  by  the  faetor 
p,  one  obtains  the  real  and  imaginary  parts  of  the  eigenvalue  7  of  smallest  real  part, 
which  must  then  be  transformed  back  into  the  original  domain.  The  procedure  from 
[19]  was  automated  for  this  dissertation  to  produce  exact  decay  rates  with  which  to 
compare  the  estimates.  While  the  procedure  just  outlined  has  already  been  developed 
in  [19],  the  numerical  results  for  specific  materials  are  new  and  show  some  interesting 
trends.  We  note  that  an  alternative  analysis  of  exact  decay  rates  has  also  been  given 
recently  in  [28]. 

Exact  decay  rates  were  computed  for  each  of  the  materials  listed  in  Table  3.1  of 
Chapter  3.  The  materials  were  initially  oriented  with  their  fibers  parallel  to  the  axial 
direction,  making  a  fiber  angle  of  0®  with  the  xi-axis.  For  this  situation,  the  strip  is 
specially  orthotropic.  See  Figure  4.1  for  an  illustration  of  the  fiber  angle  6,  for  general 
9. 


Figure  4.1:  Fiber  angle  $ 


The  principle  axes  of  the  materials  were  then  rotated  through  a  fiber  angle  range 
of  0°  —  90®.  For  values  of  0®  <  5  <  90®,  the  principal  axes  of  the  material  and 
the  axes  shown  in  Figure  4.1  do  not  coincide.  For  this  situation  the  strip  exhibits 
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anisotropy  in  the  form  of  coupling  between  stretching  and  shearing  in  the  xi  —  X2 
plane.  For  0  =  90®,  the  materials  are  specially  orthotropic,  but  with  their  fibers 
parallel  to  the  transverse  direction.  At  each  of  the  intermediate  fiber  angles,  for  which 
the  strip  is  anisotropic,  the  techniques  of  [19]  were  used  to  calculate  the  exact  decay 
rate.  A  plot  of  actual  decay  rates  is  shown  below  in  Figure  4.2.  For  convenience  we 
have  chosen  the  three  materials  HSl,  HS2  and  UM  for  which  to  plot  results.  Decay 
rate  curves  for  the  other  materials  listed  in  Table  3.1  are  similar  in  shape. 


Figure  4.2:  Exact  decay  rates  vs.  fiber  angle 


It  is  of  interest  to  note  that  for  each  of  these  materials  while  6  =  0  produces 
the  smallest  decay  rate  (as  expected),  there  is  a  fiber  orientation  angle  in  the  range 
60“  —  80“  at  which  the  decay  rate  is  largest.  This  initially  unexpected  result  has  also 
been  observed  in  other  structural  mechanics  problems  for  laminates  [20].  A  similar 
non-monotonic  behavior  of  decay  rate  with  material  properties  has  also  been  observed 
by  Choi  and  Horgan  [29]  for  a  sandwich  strip.  Mathematically,  this  peak  in  the  decay 
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rate  can  be  traced  back  to  the  existence  of  a  multiple  eigenvalue  which  was  found  in 
[19].  There,  the  real  and  imaginary  parts  of  the  modified  eigenvalue  (  were  plotted 
as  functions  of  xj)  and  u>,  forming  surfaces  in  three-dimensional  coordinate  space,  and 
a  cusp  was  found  to  exist  in  the  real  surface,  corresponding  to  a  multiple  eigenvalue 
at  that  point.  This  cusp  is  responsible  for  the  two  “kinks”  in  the  decay  rate  curves 
of  Figure  4.2.  For  each  material,  as  the  fiber  angle  is  rotated,  the  real  part  of  the 
modified  eigenvalue  (  traverses  a  path  along  the  three-dimensional  surface  mentioned 
above,  crossing  over  the  cusp  region  twice.  Due  to  the  scaling  factor  p  which  yields 
the  actual  eigenvalue  7,  the  effects  of  the  cusp  are  small  the  first  time  the  cusp  is 
crossed  and  exaggerated  the  second  time.  Thus  we  observe  in  each  of  the  curves  of 
Figure  4.2  a  small  “kink”  when  the  fiber  angle  is  near  0®  and  a  large  “kink”,  or  peak, 
as  the  fiber  angle  approaches  90®. 

Physically,  this  pealc  behavior  in  the  decay  rate  curves  may  be  related  to  the  severe 
non-monotone  behavior  of  the  elastic  parameters  €2  and  €3  with  6  in  the  particular 
fiber  angle  range  (see  e.g.  Figure  4.17  below  for  HSl).  This  peak  behavior  may  also 
occur  because  as  the  fiber  angle  increases,  the  effects  of  the  matrix  on  deformation 
become  increasingly  more  important.  Thus  for  fiber  angles  near  90®  the  material 
behavior  is  dominated  by  the  matrix  in  contrast  to  angles  near  0®  where  the  fibers 
have  the  dominant  effect  on  the  decay  rate. 

We  also  observe  from  Figure  4.2  that  for  each  of  the  materials  there  is  a  fiber 
orientation  angle  (roughly  in  the  range  30®  —  50®)  at  which  the  decay  rate  coincides 
with  that  for  an  isotropic  material  {k*  =  2.106). 
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4.6  Discussion  of  Results 

In  this  section,  we  refer  to  Figures  4.3  -  4.16  as  well  as  Table  4.1  to  discuss  the  basic 
energy  estimate,  nonlinear  optimization  estimate  and  higher-order  energy  estimates 
developed  in  this  chapter.  In  contrast  to  the  estimates  of  the  previous  chapter,  the 
estimates  of  the  present  chapter  do  not  produce  explicit  formulae  for  the  decay  rate, 
but  rather  require  numerical  methods  of  evaluation.  For  convenience  we  choose  the 
three  materials  HSl,  HS2  and  UM  as  representative  materials  for  which  to  show 
results.  We  also  wish  to  thank  Dr.  M.P.  Nemeth  of  NASA  Langley  Research  Center 
for  providing  some  of  the  numerical  codes  that  were  used  in  producing  many  of  the 
figures. 

Figures  4.3  -  4.5  show  the  exact  decay  rate  together  with  the  basic  energy  estimate 
and  the  nonlinear  optimization  estimate,  for  each  of  the  three  materials  HSl,  HS2 
and  UM  respectively.  Recall  that  these  two  estimates  were  the  only  estimates  in  the 
anisotropic  case  that  yielded  a  numerical  result  (the  higher-order  energy  estimates 
were  not  implemented  numerically  for  the  anisotropic  case).  It  is  observed  for  each 
of  the  materials  that  both  the  basic  energy  estimate  and  the  nonlinear  optimization 
estimate  (lower  bounds  for  the  exact  decay  rate)  greatly  underestimate  the  exact 
decay  rate,  particularly  in  the  middle  fiber  angle  range  of  20°  —  70°.  For  the 
fiber  angles  near  0°  and  90°,  i.e.  where  the  material  is  closer  to  being  orthotropic, 
the  estimates  are  better,  though  there  is  still  room  for  improvement.  The  nonlinear 
optimization  estimate  is  seen  to  be  slightly  better  than  the  basic  energy  estimate  2is 
the  fiber  angle  increases. 
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It  is  noted  that  in  these  preceding  figures  the  estimates  correspond  to  mathemat¬ 
ical  energy  formulations  using  the  energy  functional  Ei{z).  As  mentioned  previously, 
these  estimates  did  not  change  significantly  when  the  real  physical  strain-energy  Ep{z) 
was  used.  This  is  now  seen  in  Figures  4.6  -  4.8,  where  the  basic  energy  estimate  is 
shown  using  the  two  different  energy  formulations,  again  for  the  materials  HSl,  HS2 
and  UM,  respectively.  The  exact  decay  rate  is  reproduced  for  convenience.  While  the 
Ei{z)  formulation  is  slightly  better  in  the  outer  ranges  of  the  fiber  angle  and  the  Ep{z) 
formulation  is  slightly  better  in  the  middle  range  of  the  fiber  angle,  these  differences 
are  very  small.  Similar  results  are  seen  in  Figures  4.9  -  4.11  where  comparisons  of  the 
nonlinear  optimization  estimate  with  the  two  different  energy  formulations  axe  shown 
for  HSl,  HS2  and  UM. 

As  mentioned  earlier  in  this  chapter,  the  estimates  using  Ep{z)  have  the  advantage 
that  they  apply  to  all  physical  materials,  while  the  estimates  using  Ei(z)  require  a 
stricter  positive-definite  condition  that  may  fail  to  hold  for  some  materials.  This  can 
be  seen  from  Figure  4.12  and  Figure  4.13.  Figure  4.12  shows  the  materials  from  Table 
3.1  for  which  the  energy  functional  Ei{z)  is  a  positive  definite  quadratic  form,  i.e. 
materials  for  which  ^  ~  ^  ~  ^  >  0  (see  equation  (4.9)).  This  includes  five  of  the 
seven  materials  listed  in  Table  3.1.  Figure  4.13  shows  the  remaining  two  materials  for 
which  the  positive  definite  condition  (4.9)  fails  to  hold  at  certain  fiber  angles.  Thus, 
we  observe  that  the  Ei{z)  formulations  of  this  chapter  cannot  be  applied  to  KE  or 
BE  in  the  middle  fiber  angle  region. 

Figure  4.12  also  indicates  why  the  basic  energy  estimate  and  nonlinear  optimiza¬ 
tion  estimate  shown  in  Figures  4.3  -  4.5  yield  poor  results  in  the  middle  fiber  angle 
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range.  For  most  of  the  materials  shown  in  Figure  4.12,  the  positive-definite  con¬ 
dition  on  Ei{z)  is  barely  satisfied  in  the  middle  fiber  angle  range.  This  weakness 
of  the  energy  norm  Ei(z)  in  this  range  limits  the  accuracy  of  the  estimates  of  this 
chapter.  For  the  basic  energy  estimate  given  by  (4.40)  and  (4.45),  we  see  that  k*  is 
proportional  to  y/X  where  A  is  the  smallest  eigenvalue  of  the  positive-definite  matrix 
(4.18)  arising  in  the  quadratic  form  Wi  in  Ei{z).  It  then  follows  that  the  weaker  the 
positive-definiteness  of  Ei{z),  the  smaller  the  value  for  A  and  hence  the  poorer  the 
estimated  decay  rate  k*.  Similarly,  in  the  nonlinear  optimization  estimate  given  by 
(4.70)  -  (4.74),  we  obtained  a  result  that  was  dependent  on  the  positive-definiteness 
of  Ei{z).  Specifically,  we  found  that  in  order  for  the  conditions  (4.71)  -  (4.73)  to  be 
satisfied  for  some  choice  of  the  parsimeters  a  and  yd,  the  positive-definite  condition  on 
J5i(z)  must  hold  (see  (4.67)  -  (4.69)).  It  can  further  be  observed  that  the  weaker  the 
positive-definiteness  of  the  energy  norm  E\{^z)^  the  less  freedom  there  is  in  choosing 
the  parameters  a  and  yd  to  maximize  the  decay  estimate  and  hence  the  poorer  the 
results  expected.  Thus  the  accuracy  of  both  the  basic  energy  estimate  and  the  nonlin¬ 
ear  optimization  estimate  is  very  closely  linked  to  the  degree  of  positive-definiteness 
of  the  energy  norm  involved. 

In  Figures  4.14  -  4.16,  we  show  some  results  for  the  higher-order  energy  estimates 
of  this  chapter,  i.e.  the  quadratic  and  quartic  estimates.  For  the  anisotropic  case, 
neither  the  quadratic  estimate  nor  the  quartic  estimate  yielded  results  that  were 
very  tractable  numerically.  However,  explicit  upper  bounds  were  obtained  for  both  of 
these  estimated  decay  rates  and  are  plotted  in  Figures  4.14  -  4.16  for  HSl,  HS2  and 
UM  respectively.  The  quadratic  estimate  upper  bound  comes  from  (4.138)  while  the 
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quartic  estimate  upper  bound  is  given  in  (4.168).  For  each  of  the  materials  shown, 
we  observe  that  both  the  quadratic  estimate  upper  bound  and  the  quartic  estimate 
upper  bound  are  well  below  the  exact  decay  rate  in  the  middle  fiber  angle  range.  In 
the  outer  fiber  angle  ranges,  these  two  bounds  are  much  higher  and  in  some  cases 
exceed  the  actual  decay  rate.  Since  these  are  upper  bounds  for  the  estimated  decay 
rates  (which  are,  in  turn,  lower  bounds  on  the  exact  decay  rate),  Figures  4.14  - 
4.16  show  that  these  techniques  will  yield  very  poor  results  in  the  middle  fiber  angle 
ranges.  In  the  outer  fiber  angle  ranges  there  is  the  potential  for  accurate  results  as  the 
upper  bounds  on  the  estimates  exceed  the  exact  decay  rates,  however  the  estimates 
themselves  could  lie  well  below  the  bounding  curves.  In  fact,  for  the  HSl  and  HS2 
materials  at  15*  fiber  orientation  angle,  attempts  were  made  through  trial  and  error  to 
determine  a  numerical  result  from  the  quadratic  estimate.  The  decay  rate  estimates 
so  obtained  were  inadequate,  however.  For  these  reasons,  numerical  techniques  were 
not  further  pursued.  We  also  note  from  Figures  4.14  and  4.16  that  for  HSl  and 
UM  there  are  certain  fiber  angle  ranges  for  which  the  quartic  estimate  upper  bound 
becomes  negative.  This  implies  that  for  these  ranges  the  quartic  estimate  is  invalid 
and  cannot  be  applied. 

Results  in  the  orthotropic  limits  for  both  the  quadratic  and  quartic  estimates 
were  obtained  earlier  in  this  chapter.  The  orthotropic  quadratic  estimate  required 
that  Cl  >  .63  and  hence  this  estimate  does  not  apply  to  any  of  the  materials  in 
Table  3.1.  The  orthotropic  quartic  estimate,  however,  does  apply  to  these  materials 
and  numerical  results  are  shown  in  Table  4.1.  Here,  the  exact  decay  rates  and  the 
conservation  law  estimate  of  Chapter  3  are  included  for  comparison.  It  can  be  seen 


CHAPTER  4.  THE  ANISOTROPIC  CASE 


109 


that  for  orthotropic  materials  the  quartic  estimate  is  not  quite  as  good  as  the  conser¬ 
vation  law  estimate,  particularly  for  materials  with  larger  values  of  ci.  It  is,  however, 
better  than  the  analytic  estimate  (3.51),  excepting  the  BA  and  isotropic  materials, 
eind  better  than  (3.55)  for  small  ci. 

Finally,  in  Figures  4.17  and  4.18  we  show  typical  curves  for  the  elastic  parameters 
of  the  material  as  functions  of  the  fiber  orientation  angle.  Figure  4.17  provides  plots 
for  Cl,  C2  and  ca  for  the  HSl  material,  and  is  representative  of  most  of  the  other 
materials  used  in  this  dissertation.  As  the  fiber  angle  increases,  the  elastic  parameters 
all  vary  smoothly.  We  observe  that  ci  is  symmetric  about  45°  where  it  reaches  its 
maximum  value  while  £2  and  ca  are  symmetric  with  respect  to  each  other,  and  suffer 
a  rapid  variation  near  9  =  90°  and  ^  =  0°  respectively.  In  Figure  4.18  these  elastic 
parameters  are  plotted  for  the  KE  material  and  we  observe  some  atypical  properties. 
Here,  ci  still  varies  smoothly  and  is  symmetric  about  45°,  however  C2  and  £3  show 
sharp  jumps  in  sign.  Such  extreme  behavior  helps  to  explain  the  difficulty  involved 
with  obtaining  accurate  decay  rate  estimates  for  such  materials. 
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Figure  4.4:  Exact  decay  rate  and  estimates  for  HS2 
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Figure  4.5:  Exact  decay  rate  and  estimates  for  UM 


Figure  4.6:  Basic  energy  estimate  for  HSl  with  different  energy  functionals 
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Figure  4.7:  Basic  energy  estimate  for  HS2  with  different  energy  functionals 
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Figure  4.8:  Basic  energy  estimate  for  UM  with  different  energy  functionals 
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Figure  4.9:  Nonlinear  optimization  estimate  for  HSl  with  different  energy 
functionals 
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Figure  4.10:  Nonlinear  optimization  estimate  for  HS2  with  different  energy 
functionals 
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Figure  4.11:  Nonlinear  optimization  estimate  for  UM  with  different  energy 
functionals 
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Figure  4.12:  Materials  for  which  the  positive-definite  condition  (4.9)  is  satisfied 
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Figure  4.13:  Materials  for  which  the  positive-definite  condition  (4.9)  fails 


Figure  4.14:  Upper  bounds  (4.138),  (4.168)  on  estimated  decay  rate  for  HSl 
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Figure  4.15:  Upper  bounds  (4.138),  (4.168)  on  estimated  decay  rate  for  HS2 
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Figure  4.16:  Upper  bounds  (4.138),  (4.168)  on  estimated  decay  rate  for  UM 
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Material 

(^0 

Quartic 

Estimate 

Cons.  Law 
Estimate 

Exact 
Decay 
Rate  k* 

BE 

(0.0745) 

0.409 

0.445 

0.487 

UM 

(0.0949) 

0.306 

0.328 

0.369 

KE 

(0.1030) 

0.438 

0.467 

0.534 

HSl 

(0.1362) 

0.447 

0.468 

0.567 

HS2 

(0.1578) 

0.522 

0.564 

0.708 

SE 

(0.2366) 

0.709 

0.794 

1.195 

BA 

(0.2946) 

0.853 

1.069 

1.942 

Isotropic 

(0.500) 

0.749 

1.220 

2.106 

Table  4.1:  Quartic  estimate  for  orthotropic  materials  :  t  ~  e  « 


Chapter  5 


SUMMARY  AND 
CONCLUSIONS 


Several  methods  have  been  developed  in  this  dissertation  for  deriving  decay  rate 
estimates  for  stresses  in  an  anisotropic  linearly  elastic  semi-infinite  strip  subject  to  a 
plane  stress/strain  deformation.  The  goal  of  the  present  study  has  been  to  examine 
the  effects  of  material  orthotropy  and  anisotropy  on  the  decay  rate  of  local  stresses 
generated  by  self-equilibrated  edge  loadings. 

The  problem  was  analyzed  first  for  the  important  subclass  of  specially  orthotropic 
materials,  where  two  estimates  (the  analytic  estimate  and  the  conservation  law  esti¬ 
mate)  were  obtained  that  provide  lower  bounds  on  the  decay  rate  for  stresses  in  the 
material.  The  analytic  estimate  gave  an  explicit  formula  in  terms  of  nondimensional 
material  parameters  and  thus  directly  revealed  the  effects  of  material  orthotropy  on 
the  stress  decay  rate.  Two  essential  parameters  required  for  determining  the  stress 
decay  rates  have  been  identified  as  the  dimensionless  parameter  Ci  and  the  beta  ratio 
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(^)i.  Upon  expressing  these  parameters  in  terms  of  the  usual  engineering  material 
constants,  the  stress  decay  behavior  can  be  found  2ind  trends  can  be  identified  for  a 
wide  class  of  materials.  It  has  been  shown  that,  in  general,  materials  with  the  smallest 
values  of  ei  have  the  smallest  decay  rates  and  hence  the  largest  decay  lengths.  When 
compared  with  the  corresponding  exact  stress  decay  rates,  the  analytic  estimate  has 
been  shown  to  preserve  relative  ordering  of  the  decay  rates  among  the  various  mate¬ 
rials  from  Table  3.1  and,  although  an  underestimate  for  the  exact  decay  rate  (and  a 
overestimate  for  the  exact  decay  length),  the  analytic  estimate  yields  a  conservative 
and  useful  result  that  captures  the  overall  behavioral  trends  of  specially  orthotropic 
materials.  While  the  conservation  law  estimate  did  not  produce  an  explicit  formula, 
numerical  results  were  obtained  that  axe  significantly  better  than  the  corresponding 
results  obtained  using  the  analytic  estimate,  and  for  some  materials  were  in  good 
agreement  with  the  exact  stress  decay  rates.  The  improved  accuracy  of  the  conser¬ 
vation  law  approach,  particularly  for  the  strongly  orthotropic  materials  considered, 
suggests  the  usefulness  of  such  a  method  for  estimating  solutions  to  similar  problems. 
For  both  the  analytic  and  conservation  law  estimates,  an  asymptotic  analysis  resulted 
in  asymptotic  formulae  for  the  estimated  decay  rates  for  materials  with  small  values 
of  €i.  For  strongly  orthotropic  materials  these  asymptotic  formulae  resulted  in  stress 
decay  rate  estimates  that  are  proportional  to  a  result  that  has  already  been 

applied  to  many  practical  problems  involving  composite  materials.  It  hais  been  fur¬ 
ther  shown  that  for  all  of  the  analyses  presented  here  for  the  orthotropic  strip,  the 
most  accurate  decay  rate  estimates  are  obtained  for  materials  with  small  values  of 
the  orthotropic  parameter  Ci. 
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The  strip  problem  has  also  been  examined  for  anisotropic  materials,  where  the 
difficulty  encountered  in  obtaining  results  greatly  increases.  For  this  problem,  four 
decay  rate  estimates  have  been  derived;  i.e.,  a  basic  energy  estimate,  a  nonlinear 
optimization  estimate  and  two  higher-order  energy  estimates.  Unlike  the  analyses  for 
the  orthotropic  strip,  explicit  decay  rate  formulae  for  the  anisotropic  strip  in  terms 
of  the  elastic  material  constants  were  not  obtained  due  to  the  increased  analytical 
complexity  caused  by  the  presence  of  anisotropy.  Both  the  basic  energy  estimate  and 
the  nonlinear  optimization  estimate  yield  numerical  results  that  depend  implicitly 
on  the  elastic  parameters  of  the  material,  ei,  €2,  €3  in  addition  to  the  beta  ratio. 
These  estimates  have  been  compared  with  corresponding  exact  decay  rates  computed 
numerically  for  a  set  of  specially  orthotropic  materials  in  which  the  principal  material 
axes  were  rotated  through  a  fiber  angle  varying  from  0°  to  90°.  Both  the  basic  energy 
estimate  and  the  nonlinear  optimization  estimate  exhibited  similar  characteristics; 
they  gave  more  accurate  estimates  for  the  fiber  angles  near  0°  and  90°,  where  the 
material  is  close  to  being  orthotropic,  while  giving  less  accurate  results  in  the  middle 
fiber  angle  region  where  more  general  anisotropy  is  observed.  The  limitations  of  both 
of  these  approximate  methods  for  anisotropic  strips  has  been  found  to  be  directly 
related  to  the  positive-definiteness  conditions  on  the  energy  norms.  The  higher-order 
energy  estimates  were  modelled  after  the  conservation  law  approach  of  Chapter  3 
2ind  resulted  in  optimization  problems  which,  due  to  the  number  of  undetermined 
parameters  involved,  did  not  result  in  tractable  numerical  results.  These  higher-order 
estimates  were  found  to  simplify  considerably  in  the  orthotropic  limit.  However  the 
numerical  results  obtained  were  not  better  than  the  orthotropic  estimates  of  Chapter 
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3.  Asymptotic  analyses  were  not  pursued  for  the  anisotropic  strip  since  small  values  of 
the  parameters  C2  and  €3  do  not  have  the  physical  significance  that  small  values  of  61  do 
for  the  orthotropic  strip,  i.e.  corresponding  to  strong  orthotropy.  In  particular,  small 
values  of  62  and  €3  are  not  exhibited  by  common  materials  of  technological  interest. 
In  general,  the  methods  of  this  dissertation  when  extended  to  the  anisotropic  strip 
were  found  to  be  difficult  to  implement  and  were  not  able  to  provide  as  accurate 
results  as  desired.  Nevertheless,  the  results  indicate  some  of  the  difficulties  involved 
in  modelling  the  anisotropic  strip  and  give  insight  into  what  may  be  done  to  improve 
the  estimates.  Furthermore,  since  the  results  yield  lower  bounds  for  the  decay  rates 
and  thus  upper  bounds  for  the  decay  lengths,  the  results  are  conservative  and  thus 
immediately  useful  from  an  engineering  design  perspective. 

We  conclude  with  an  outline  of  some  area.s  of  future  research.  Ideally,  an  explicit 
decay  rate  formula  in  terms  of  €1,  €2  and  €3  is  desirable  for  the  anisotropic  strip.  This 
would  clearly  reveal  the  importance  of  anisotropy  on  the  decay  rate  of  edge  effects 
and  would  be  applicable  to  a  broad  range  of  materials.  This  insight  could  stimulate 
the  investigations  of  many  other  problems  using  a  classical  approach  with  anisotropic 
elasticity.  Such  insight  could  greatly  eissist  in  tailoring  composite  structures  and 
provide  guidance  for  modelling  complex  structures  with  powerful  computational  tools 
such  as  the  finite  element  method.  One  promising  approach  to  solving  the  anisotropic 
strip  problem  is  the  conservation  law  method.  While  one  of  the  two  conservation  laws 
developed  in  this  research  was  shown  to  yield  highly  accurate  results  for  orthotropic 
materials,  it  remains  to  be  seen  whether  the  remaining  conservation  law  can  be  used 
in  a  similar  fashion  or  if  additional  conservation  laws  can  be  developed  to  improve  the 
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results.  Finally,  as  the  anisotropic  analysis  revealed,  the  decay  rate  estimates  obtained 
here  are  not  accurate  whenever  the  positive-definiteness  of  the  energy  norm  is  weakly 
satisfied.  This  suggests  that  alternate  energy  norms  might  be  more  appropriate  for 
the  analysis  of  the  anisotropic  strip. 
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